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1.0 INTRODUCTION; 

This report summarizes progress made during the past year and 
a half on a research program funded by NASA/LaRC and conducted by 
Prof. D. Killinger's Lidar Research Group at the University of South 
Florida. This research has been directed toward the development of 
eye-safe, solid-state Lidar systems, with an emphasis on Coherent 
Doppler Lidar for Atmospheric Wind Measurements. Details of the 
Tasks and progress made are in the following sections of this report. 
The most important details of the progress made are presented in the 
following sections, along with supporting related documentation 
(reprints of scientific journal articles) which are enclosed in the 
Appendix of this report. 

1.1 OVERVIEW OF PR OGRAM TASKS: 

The original program had Tasks to develop the Ho Lidar system, 
compare its performance with a C02 Lidar system for the detection 
of atmospheric backscatter, and develop (with NASA/LaRC) a 
coherent Doppler Ho Lidar system for potential windshear detection. 
This program was divided into two phases. Phase I (Task 1) 
involved the measurement of atmospheric backscatter coefficient, [J, 
for the Ho and C02 Lidar systems. During this time, it was 
anticipated that NASA/LaRC would be finishing the development of a 
single-frequency Ho laser source. Phase II (Task 2) involves the 


cooperative development of a coherent Doppler Ho Lidar using the Ho 
laser developed at LaRC and the Ho Lidar system at USF . 

Task 1, involving the measurement of the backscatter coefficient 
of the atmosphere, p , has been completed, and Task 2, involving the 
development of a Ho Coherent Lidar, has been initiated. A slight 
modification to the original program was made to better 
accommodate the schedual of the completion of the single-frequency 
Ho laser and to allow involvement in NASA/Lockheed atmospheric 
wind measurements during a recent Shuttle Launch at KSC. As such, 
an additional subtask involving measurements of atmospheric 
refractive turbulence was made in support of another NASA Lidar 
program, instead of the continued CO2 Lidar studies. 


2.0 SPECIFIC PROGRESS ITEMS; 

Significant progress has been made already under this program. 
We have recently demonstrated for the first time (1) the tuning of 
the Ho laser over a water vapor absorption line in the atmosphere , 

(2) measured simultaneously the range-resolved Lidar mapping of 
atmospheric aerosols and water vapor concentration, (3) measured 
the temporal and range-resolved attenuation and backscatter 
coefficient of the atmosphere, (4) measured the temporal variation 
in atmospheric refractive turbulence and performed related Coherent 
Lidar performance calculations, and (5) made initial progress in the 
joint development of a single-frequency, high power Ho 
Laser/Coherent Lidar system. Details of these items are presented 
in the following sections. 

2.1 Tunable Ho DIAULidar Atmospheric Measurements 

The direct-detection 2.1 pm Ho:YSGG Lidar system was developed 
and used to measure the range-resolved backscatter coefficient of 
the atmosphere. In addition, through the use of intracavity tuning 
etalons, the wavelength of the Ho laser was tuned across an 
absorption line of water vapor in the atmosphere; range-resolved 


Lidar measurements of atmospheric water vapor were also 
demonstrated. 

Figure 1 shows initial Lidar/DIAL Ho measurements of the 
atmospheric backscatter and water vapor density. As can be seen, 
even with a low power 7 mJ/pulse Ho laser, range resolved 
measurements at ranges of several kilometers is easily obtainable. 
Details of the development of the Q-switched Ho laser source are 
given in the Appendix (reprint of Journal article), and details 
concerning the initial water vapor measurements of the atmosphere 
are also given in the Appendix. 

An improved version of the Ho Lidar system has recently been 
completed and enhanced performance of the system has also been 
demonstrated. In thios case, the Ho laser was upgraded to include 
better output coupling mirrors, an A/O Q-switch, and low loss 
transmitting mirrors. The output parameters of the Ho laser was 25 
mJ/pulse and a 90 ns pulselength. Figure 2 shows a schematic of 
the improved Ho direct-detection Lidar system. Range-resolved 
measurements have been obtained out to a range of 3 km. 
Measurements of the temporal variation of the atmospheric 
attenuation coefficient, a, and the backscatter coefficient, p, were 
made, inorder to calibrate the system sensitivity and to better 
understand the important Lidar atmospheric parameters for a 2.1 pm 
Ho Lidar. Figure 3 shows our Lidar data obtained over a 24 hour 
period. The "error bars" do not represent the accuracy of the 
measurement, but show the range-resolved variability in the value of 
a or p measured along the line-of-sight of the Lidar beam; the range 

resolution of the system was about 100m. 

Of interest is the data shown in Fig. 4, which presents the ratio of 
a/p. This ratio shows a large change during the day and seems to 
indicate that the composition of the backscatter aerosols is changing 
from predominately water based aerosols to dry-particulate based 
aerosols. The ratio of a/p on the order of 10 to 100 tends to indicate 

water vapor aerosol formation , and a ratio on the order of 100 to 
500 tends to indicate dry or dust aerosol composition. Further 
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Fig. 1 Ho lidar/DIAL measurements of the temporal vnrint ion or (a) the nfmospherir aerosol backscaUering profile and (b) (hr almospherir 
water -vapor concentration profile. 'Fhe lidar was pointed nt approximately a 4-2" elevation angle. 
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Fig. 3 Diurnal and Range-resolved variation in Ho Lidar 
measurements of atmospheric extinction coefficient and 
backscatter coefficient. 






measurements are being planned to better quantify these results; a 
journal paper has been submitted on these results. 


2.2 Atmospheric Turbulence Measurem ents and D etector Array s 

Early work which we performed in 1989 showed experimentally for 
the first time that atmospheric refractive turbulence would severely 
hurt the performance of some Coherent Lidar systems operating in 
the near to mid-IR; this work was part of another study which 
investigated a 1 pm Coherent Doppler Lidar for rocket plume 
measurements. The results of these studies were applied to the case 
of a 2.1 pm Ho Lidar for atmospheric wind sensing. The pertinent 
results are shown in Fig. 5 which gives the predicted maximum Lidar 
telescope size as a function of wavelength, range, slant path angle, 
and turbulence level; here WN (Winter Night) refers to a turbulence 
level of C n 2 - 2 x 10 -12 nrr 2/3 , and SD(Summer Day) to C n 2 - 1 * 
10-13 m -2/3. As can be seen, a Ho Coherent Lidar may be limited 
to telescope apertures on the order of 20 - 50 cm under high 
turbulence levels and horizontal (ground based) operation. Figure 6 
shows similar data for a 2 pm Ho Coherent Lidar operating from 
space. As can be seen, rather large telescopes are usable. Details 
of these studies are presented in the Appendix and in another journal 
article which has been submitted for publication. 

It should be noted that the reduced Coherent Lidar signal due to 
atmospheric refractive turbulence can be overcome by use of a 
multi-element heterodyne detector array. We have recently shown 
that the Lidar S/N can be improved by a factor of N 2 through use of a 
N x N heterodyne detector array. ^ We plan to investigate the use of 

such an array for the enhancement of the Ho Coherent Lidar S/N, and 
the use of the array for target tracking and "Lag Angle" compensation. 
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Fig. 5 Predicted maximum useful Lidar telescope size for a Coherent 
Doppler Lidar as a function of wavelength, slant angle, Lidar range, 
and turbulence level. 
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Fig. 6 Predicted maximum Coherent Lidar telescope size for a 2.1 
pm Ho Lidar in a satellite orbit and looking downward. 
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2.3 Diurnal Measurements of Cp ^ - for KSC Udar Mf&smgJUSOl S 

A C n 2 scintillometer instrument was obtained on loan from NOAA 
and used to measure the level of atmospheric refractive turbulence. 
This instrument was used to compare and validate our 1 pm Coherent 
Lidar measurements and the Ho direct-detection Lidar measurements. 
Recently, we used this instrument in support of a related 1 pm 
Coherent Doppler Lidar Wind-Shear measurment test by NASA/LaRC 
(D. Moerder) and Lockheed ( R. Targ and J. Hawley) at KSC during a 
recent (Sept. 12, 1991) Shuttle Launch. The values of Cn 2 obtained 
are shown in Fig. 7. As can be seen, the value of the turbulence level 
varies dramatically during the day. Fig. 8 shows the predicted 
maximum useful telescope aperture for the upward-looking 1 pm 
Lockheed Lidar system using the data of Fig. 8. As can be seen, the 
20 cm diameter Lockheed Lidar was probably not affected by 
atmospheric turbulence, due to its upward (verticle) looking Lidar 
beam direction. Details of these measurements are contained in a 
Masters Thesis^ and are being prepared for publication. 


2.4 Joint Development of Sing le- Frequency Ho Laser/Lida r 

A joint program has been initiated with NASA/LaRC (M. Storm, N. 
Barnes) and USF in the construction/development of a single- 
frequency Ho laser/Lidar. NASA has developed a preliminary version 
of the laser and delivered it (on loan) to USF for system testing and 
laser modifications. The Ho laser consists of a diode-pumped Ho.YAG 
CW laser and a Q-switched, flashlamp pumped Ho:YAG amplifier. We 
have recently been in the process of characterizing this laser and 
modifying the cavity design, in consultation with M. Storm. Our 
preliminary results indicate that the diode-pumped Ho:YAG laser 
operates multi-frequency on the 2.09 pm Ho line, but operates 
single-frequency on the 2.12 pm line; this is done by tuning the 
temperature of the Ho crystal. We are currently working on 
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Fig. 8a Predicted Diurnal Variation of Coherence Length, p 0 . 
and Effective Telescope Diameter, D u , lor a 1pm 
Coherent LIDAfl (Data for 9/12/91). 



Fig. 8b Measurements Showing the Diurnal Variation of Cn 2 
at KSC for 9/12/91. 
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Fig. 9 Output power of modified NASA/LaRC Ho:YAG flashlamp laser 
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modifying the diode-pumped Ho laser to operate single-frequency on 
the 2.09 pm line, since this line has higher gain than the 2.12 pm line. 

The flashlamp cavity has been considerably modified by using a 
different optical mounts and cavity assemply. The output of the 
flashlamp Ho laser is shown in Fig. 9, which indicates that 
approximately 24 mJ/pulse, TEMoo> Q-switched laser pulses were 
obtained; close to 100 mJ/pulse was obtained for non TEMoo 
operation. Work is now being conducted in injection seeding the 
flashlamp Ho:YAG laser using the 2.12 mm line; operation at 2.09 mm 
will also be studied. Future plans include using this laser in our Ho 
Lidar system and the operation of the system for Coherent Lidar 
measurements of the atmosphere. It is anticipated that Lidar 
measurements of atmospheric winds will be made inorder to quantify 
the potential sensitivity of the Ho Doppler Lidar for wind detection. 


3 ANTICIPATED FUTURE PROGRAM TASKS 

Future program tasks for FY91-FY92 will include (1) continued 
development of single-frequency Ho:YAG laser source in conjunction 
with NASA/LaRC, (2) incorporation of Ho:YAG laser into a Coherent 
Doppler Lidar system, (3) investigation of use of 2-D heterodyne 
detector arrays to improve Ho Coherent Lidar performance, and (4) 
joint development/incorporation of Ho:YLF laser source into Ho Lidar 
system. This latter task is related to the current NASA/LaRC 
development of a high-power, diode-pumped Ho:YLF laser source for 
atmospheric Lidar. 

It is anticipated that the Ho.YLF laser may have to be slightly tuned 
in wavelength inorder to increase its transmission through the 
atmosphere. This can be seen in Fig. 10, which is a plot of the 
transmission of the atmosphere over a 10 km path (ground level) and 
the overlap of the tuning range of the lines for Ho:YAG, Ho:YSGG, and 
Ho:YLF. Futher studies will be made inorder to verify the above 
predictions with actual DIAL/Lidar atmospheric measurements. 
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Fig. 10 Atmospheric transmission through 10 km path (sea level) 
showing overlap of tuning range of Ho:YAG, HoiYSGG, and Ho:YLF laser 
lines. 
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Appendix 1 (Reprint of article in Japanese Journal of Optics, 20,612 (1991) 


Performance Characteristics of Acousto-Optic Q-switched 
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2. Experimental Setup 

The laser used in this experiment was a com- 
mercially available, lloshlamp pumped. Ho: YSGG 
laser (Schwartz Electro-Optics, baser 1 2-3) ; a srhe- 
malic diagram of the laser is shown in 1 iff* 1* 1 he 
laser cavity used a 3 min X 75 mm Ho: YSOG rod 
with Cr and Tin co-doping to aid in Hashlamp 
operation; the concentration was fix 10 n cm 1 Ho, 

2 x 10' 1 cm" 1 Cr, and 8: 10"° cm' 3 Tin. The laser 
cavity had a 100% rdicci ion rear mirror, 00% re- 
flection output mirror, and a cavity length of 07 
cm. A fused quartz A/O Q-swilch (Newport Elec- 
tro-Optics Systems, Model No. N .12027 .;0 
2 . 1) was used in the laser cavity and was modu- 
lated with a radio frequency (RE) signal of 27. 1 
MHz supplied by an EE driver (Newport Electro- 
Optics Systems, Model No. N 31027 50 DM). The 
A/O device was placed between die rear mirror and 
the laser rod. An aperture was inserted between 
the A/O device and the rod to achieve TKM«n single 
spatial mode laser operation. I be pulse encigv 
was measured with an energy meter and the shape 
of the Q-swilch pulse was monitored using a fast 
rise-lime InGaAs (Epilaxx, Model No. ETX 1000 
OR 21) pholodctcclor with a high speed amplifier 
(Analog Module, Model No. 313 1 S). The wave- 
length of the laser was monitored by a 0. 3 m grat- 
ing spectrometer. 

The relative timing of the Hashlamp, the build- 
up of the population inversion nf the 1 In atoms, and 
the Q-swilch trigger pulse arc important, in deter- 
mining the output characteristics of the Ho lasci 
In order to measure some of these temporal param- 
eters, a Go detector (with a Gc window) was used 
in the cavity in order to measure the 1.8 //in fluo- 
rescence from the llo laser rod, a Si detector was 
used adjacent to the Hashlamp to measure the visi- 
ble output of- the Hashlamp, and another Si de- 
tector was used to measure the green fluorescence 
from the laser rod. An InGaAs detector monitored 


the laser oulpul pulse. 

3. Results 

Figure 2t a ) shows a composite photograph nf 
the relative timing of the detected optical and Her- 


f a i 


( I. 


Rig. 2 fa) Photograph "f die relative 'mum: 
nf the Hashlamp trigger. Hashlamp visible output. 
Ho fluorescence near 1. 8/uu, Hn fluorescence in 
(he green wavelength range, normal mode osol- 
j n t ion of (he llo laser (without A/O Q-s witch), 
RR f |,ive signal to the A/O Q-swilch, and «hr 
Q-sw itched laser pulse, (h) Photograph of the 
pulse sliape of the Q-switrhed lh» laser. E.yrr 
conditions arc Hashlamp input energy of I2IJ. 
f Jo laser mil temperature of 13 th and repehlmn 

rate 1 1 f d I f z . 
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U i r; ,i signals. '['lie signals show, in order from top 
to bottom, tlie relative liming of tlie llasldamp trig- 
ger, llasldamp visible output, Ho fluorescence new 
| 8 /mu, Ho fluorescence in the green wavelength 
range, normal mode oscillation of tbe Ik* Insci 
( without U-sw ilclO, tbe l\K drive signal to tbe A/O 
Q-switch, and tlie (^-switched 1 fo laser pulse, ll 
should be noted dint tbe signal obtaiue<l from the 
|.8 /mu No fluorescence is indicative of tbe mea- 
sured T bins lifetime of ibis level and is similar to 
1 1 ml repotted previously/’* 7 ’ 1 be strong green flu- 
orescence is due to the r, S 2 - : ’k transition in Ho/’ 
Figure -(b) shows a more detailed pulse shape 
of the O- switched No laser pulse. 

The pulse energy and pulse width of tbe Q- 
sw itched laser pulse were found to be strongly de- 
pendent upon tbe delay time of the A/O Q switch 
RF i urn-olf time. 'This is shown in Fig. d which 
shows the measured pulse energy as a function of 
the delay time. 'Hie optimal delay lime was found 
to be approx i m ally G20 /rs for our experimental 
setup- Our laser had the following operational 
parameters : l KM no mode, 12 mj laser pulse en- 

ergy, 120 ns pulse length, llasldamp input energy 
of 121 J; pulse, and No rod temperature of 10'C* 
An aperture of 2 mm was used to force the laser 
output into a TKMiv mode. The pulse energy for 
O- switched mullimode operation was 20 mj/pulse 
and for normal oscillation (no Q-swilch) was 27 
mj/pul.-e. Laser threshold was 9G J for the flash- 
lamp energy. As seen in big. 2, the laser output 
is maximum al a pulse repetition frequency (I’lvh) 
uf 2 II/. with a decrease seen al 1 1 h ; a significant 
decrease was also observed at b Nz due to an in- 
crease hi the thermal lens ellecl al the higher KMT 
values. 

The (0- switched pulse width was also measured 
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as a fundkTn^of the cT* witch delay lime and is 
shown \f Fig. I. As/‘ an he seen, the shortest 

pulse Itun pit occm '^fffuie delay time for the high- 
est output power which is as expected. Il should 
be added that because the transfer of excitation in 
the laser rod involves the Cr and Tin ions also, the 

1 lo -population inversion can also occur at signifi- 
cant delay lime. 10 This long term excitation pro* 
eess can produce a secondary (or double) pulse 
which occurs much later in time aller the Q-switch- 
t-d pulse. In order to clliminale tins double pulse, 
the kl* drive l<> the O switch was turned "back- 
on” after apjiroximalely Lb/rs at the end of the 
regular RF drive pulse. 'Phis is not shown in Fig. 

2 for ease of clearity in explaining the relative 
timing of the* No laser. 

The output energy and pulse width of the No 
laser were also investigated as a function of the 
llasldamp input power and operating temperature 
of the I lo laser rod. Figure b shows tlie measured 
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laser output pulse energy for both normal oscilla- 
tion anti Q-switchcd operation, and the measurer 
Q-switehed pulse width. As seen in Fig. 5, at the 
higher llashlamp input levels, the TEMoo mode Q- 
sw itched laser pulse had a pulse energy of 25 nij 
and a pulse width of 88 ns with a Hash lamp input 
energy of 156 J. The normal oscillation pulse en- 
ergy increased to nearly 150 mj/pulse at the higher 
input levels but the Q-switchcd output power did 
not significantly increase. While this may be in- 
dicative of the relatively low saturation llucncc of 
the Ho loser, it may be also due to losses in our 
A/O Q-switch. We arc currently studying this 
more in ordpf^tSStnprovc the Q-swilehed output 
power. yS J 

Figure 6 sh£<vs the Q-switchcd pulse energy 
and/pulse wktfhas a function of the laser rod tern- 
perahw^As can be seen, an operational tempera- 
ture near 16"C seems to optimize the output power 
of the Ho laser with a concurrent reduction in the 
laser pulse width. It may be added that, in gen- 
eral the threshold energy of the Ho laser is decreas- 
ed when the rod is cooled because the population 
of the lower state laser level of Ho is lowered. 9 
The optimal temperature is dependent upon several 
experimental conditions, including the rod geome- 
try, cooling design, related thermal aspects, heat 
load of the llashlamp, and laser energy extraction 
efficiency. In our experiments, the optimal tem- 
perature was found to be near 16T. although this 
value will change slightly for different PKF and 
tlnshlamp input energy. 

We have also studied the inability or the 1 lo 
laser. Two uncoatcd glass clnlons (0. 2 mm and 1 . 0 
mm thickness) were inserted in the Ilo laser cavity 
between the output mirror and the laser rod. I lie 



Fig. 6 Q-switched pulse energy and Q-switched 
pulse width as a function of the laser rod 
temperature. 
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measured tuning curve is shown in 1 tg- 7 which 
indicates that the Ho laser can be tuned smooth y 
from 2. 081 /tm to 2. 087 /mi ; the linewidth of the 
laser was measured to he approximately 0. 1 nm. 
These tuning results arc similar to our previously 
reported results." It is informative to mention 
that the wavelength of the Ho : YSGG laser is s.m- 
tlar, but not exactly the same, lo that of the Ho : 
YAG laser. Recently, Henderson and Hale"” have 
shown that the Ho: YAC. laser operates on two 
closely spaced lines at 2. 090 /mu and at 2. 097 /mu, 
with tuning obtained from 2.088 to 2. 098 /mu. 


net/1 n 


In this technical note, we have documented our 
studies of the operational characteristics of a flash- 
lamp pumped A/O Q-switchcd Ci . Tm . He 
YSGG laser operating on the 2. 1 /mi I lo laser line 
Typical operating parameters of the Q switched 
laser were pulse energy of 25 m. I and a pulse width 
of 88 ns with a flashlnmp input energy of MiJ. 
temperature of IfTC. ami repetition rale of .Hlz. 
The pulse energy was somewhat greater t inn l ml 
obtained previously using a spinning mirror Q- 
switch" Continuous tuning of the Ho laser was 
also obtained over a range from 2. 081 /mu to 2. 087 


/nn. . .it 

Our future plan is to injection seed our Ho 

loser in order lo have .1 operate as a tunable single- 

frequency laser source. Such operation will permit 

it to be used in a Doppler lidar system as well as 

a differential-absorption lidar Tor atmospheric water 

vapor and CO 2 measurements. 
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Tunable 2.1-^m Ho lidar for simultaneous range-resolved 
measurements of atmospheric water vapor and aerosol 
backscatter profiles 


Sungdo Cha, Kin Pui Chan, and Dennis K. Killinger 


An eye-safe, tunable differential-absorption lidar system has been developed for the range-resolved measure- 
ment of aerosol backscatter and water vapor in the atmosphere. The lidar uses a flash-lamp-pumped, Q- 
switched, 10-mJ solid-state Ho:YSGG laser that is continuously tunable over a 20-cm' 1 wavelength range 
near 2.084 pm. Both path-averaged and range -resolved measurements were performed with the Ho differen- 
tial-absorption lidar system. Preliminary measurements have been made of the temporal variation of 
atmospheric aerosol backscatter and water-vapor profiles at ranges out to 1 km. These results indicate that 
the Ho lidar has the potential for the eye-safe remote sensing of atmospheric water vapor and backscatter 
profiles at longer ranges if suitably enhanced in laser power and laser linewidth. 


I. Introduction 

Recently there has been increased interest in the de- 
velopment of an eye-safe, solid-state laser-lidar sys- 
tem for the detection of wind shears in front of aircraft 
and for the remote sensing of water vapor in the atmo- 
sphere. 1 3 The 2.1-jum Ho laser is considered to be 
potentially well suited for these lidar applications be- 
cause it is eye safe (X > 1.4 jxm), can be pumped by 
diode lasers, and offers the potential for tunability/ 1 
Previously we reported the development of a Ho 
lidar and the range-resolved remote sensing of atmo- 
spheric aerosol backscatter profiles at ranges out to 2 
km. r> In this paper we report an extension of this work 
that includes the precision tuning of our low-power (10 
mJ/pulse) Ho laser over a 20-cm" 1 wavelength range. 
The continuous-tuning coverage was measured from 
2.080 to 2.089 ^m (—90 A) and also from 2.104 to 2.106 
M m. The wavelength range of 2.080-2.089 overlaps 
several absorption lines of water vapor so that the Ho 
laser can be used for the differential-absorption lidar 
(DIAL) remote sensing of water vapor in the atmo- 
sphere. As an initial demonstration of the potential 
utility of the Ho lidar-DIAL system, the atmospheric 
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aerosol backscatter and water-vapor concentration 
profiles were measured simultaneously at ranges out to 
1 km as a function of time during the day. This is to 
our knowledge the first DIAL measurement of atmo- 
spheric water vapor made with the Ho laser and shows 
that the Ho DIAL system may be potentially well 
suited for the detection of both aerosol backscatter and 
water-vapor concentration profiles in the atmosphere. 
Our experimental results indicate that future exten- 
sion of the Ho laser to higher power, of the order of 1 J/ 
pulse, should provide up to tens of kilometers detec- 
tion range. 

It should be added that our range-resolved lidar 
measurements of atmospheric water vapor and aerosol 
backscatter represent an extension of previous work by 
us and by others. Previously Baker, 1 R Hardesty, 7 and 
Grant et a/. 8 measured water-vapor profiles by using 
CO 2 lasers, while Browell et a/. 9 and Ehret et a/. 10 used 
dye lasers tuned in the near-IR spectral region. Of 
course, early research by Schotland 11 and by Zuev et 
a/. 12 showed the first use of a solid-state (ruby) laser 
source for DIAL measurements of water vapor, as did 
our recent research that used a Co:MgF 2 solid-state 
laser. 13 What is novel about our Ho DIAL research 
presented in this paper is that the Ho laser is believed 
to offer several advantages over these previous laser 
sources, namely, that it can be pumped by diode lasers, 
it operates in the eye-safe wavelength range, and the 
atmospheric backscatter coefficient at 2.1 ^m is great- 
er than that at 10 pm. Thus the Ho lidar has the 
potential to be important not only for atmospheric 
aerosols and water-vapor profiling but also for three- 
dimensional Doppler mapping of wind velocity. 
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II. HoiYSGG DIAL System 

In this section we describe wavelength tuning of the Ho 
laser by an intracavity etalon and the Ho DIAL system. 
The Ho DIAL system consisted of a tunable Ho:YSGG 
laser, spectrometers to monitor the laser wavelength, 
lidar optics (a beam-steering mirror and a receiving 
telescope), optical detectors, and a data-acquisition 
system to record the lidar return signals. A schematic 
of the Ho DIAL system is shown in Fig. 1. This Ho 
lidar is essentially the same as reported previously 5 
with the addition of the laser-tuning elements and 
associated wavelength monitors. 

A. Wavelength Tuning of the Ho Laser 

Previous research on the 2.1-pm Ho laser has shown 
that the Ho laser can be used for tunable DIAL appli- 
cations. First the Ho laser has been shown by Johnson 
et al. to operate on several different emission lines over 
a wavelength range of —2.09-2.12 pm. 14 Barnes and 
Gettemy showed that a grating could be used to tune 
the Ho:YAG laser over each of these emission lines 
with a tuning range of —30 cm” 1 (i.e., 0.013 pm). 15 
More recently Henderson and Hale have shown diode 
laser-pumped, single-frequency operation of a cw 
Tm:Ho:YAG oscillator, which indicated the potential 
of the oscillator as a narrow spectral laser source. 16 
Our development of the tunable Ho laser used a com- 
mercial flash-lamp-pumped Ho laser (Schwartz Elec- 
tro-Optics, Model Laser 1-2-3) with a Ho:YSGG laser 
rod; HorYSGG was used instead of Ho:YAG because of 
its immediate availability. 

When the laser was operated in a non-Q-switched 
mode, its output pulse length was 250 ps with 20-mJ 
pulse energy. In the Q-switch operation using a spin- 
ning mirror, the pulse length was 200 ns with — 10-mJ 
pulse energy. The pulse repetition rate was limited 
primarily by the spinning mirror Q switch and was 2 
Hz. The wavelength of the laser output was moni- 
tored by a 0.5-m spectrometer and a 2.1 -pm sensitive 
TV video camera, which is described in more detail 
below. 

In the initial experimental setup, the tuning of the 
Ho laser was accomplished by using a single, uncoated 
glass §talon in the laser cavity. The laser wavelength 
was changed by varying the tilt angle of the etalon. 
The measured tunable wavelength range of the 
HorYSGG laser using the intracavity etalon (thickness 
0.17 mm) is shown in Fig. 2(a). The wavelength of the 
Ho laser was tuned continuously from 2.080 to 2.089 
/im in single-transverse-mode (TEMoo) operation and 
was also tuned near the 2.105-pm line in multimode 
operation. The latter line has a lower gain than the 
2.08-pm line but shows greater output power in Fig. 
2(a) because multimode operation uses a much larger 
active mode volume in the laser rod than does the 
TEMoo mode. The tuning range of the multimode Ho 
laser was measured from 2.1040 to 2.1064 pm. This 
tunable spectral range is consistent with previous fluo- 
rescence emission data for Ho. 17 It should be noted 
that single, uncoated glass ^talons with thicknesses of 
0.1 and 0.2 mm, respectively, were also tried; the wave- 



Fig. 1. Schematic of the Ho lidar-DIAL system; Det, detector; 
pyro, pyroelectric. 
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Fig. 2. Measured tuning range of the Ho:YSGG laser using tilted 
intracavity Stalons. (a) One talon (thickness 0.17 mm) was used in 
the non-Q-switched laser cavity, (b) Two Stalons (thicknesses 0.2 
and i.O mm) were used in the Q-switched laser cavity. 
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length tuning range was measured to be —80 and 70 A. 
By using the 0.17-mm etalon, the linewidth of the laser 
was reduced from 3 to 0.46 cm" 1 . 

In our initial experiments it was noticed that the 
tuning of the Ho-laser wavelength was not completely 
smooth but occurred in a series of etalon mode hops 
with a spacing between wavelengths of —3 A (0.68 
cm" 1 ). It was determined that the output coupler of 
the Ho laser was acting as an Etalon. The output 
coupler was made of ZnSe (thickness 3 mm, reflectivity 
70%) with flat, plane-parallel surfaces and had a calcu- 
lated free-spectral range as an etalon of —3 A. Several 
different thicknesses and types of output coupler were 
also tried, which confirmed this hypothesis. As a re- 
sult, a 1° wedged output coupler was used and pro- 
duced a smooth tuning of the laser wavelength within 
the accuracy of our spectral-measurement resolution. 

To reduce the linewidth of the Q-switched Ho laser 
even further, two uncoated glass etalons (thicknesses 
0.2 mm and 1.0 mm) were mounted on a single rotating 
stage and placed in the Ho-laser cavity. To eliminate 
mode hopping over the transmission modes of the 
thicker Etalon, the transmission peaks of the two 6ta- 
lons were adjusted to the same position in the wave- 
length domain. This was done by careful adjustment 
of the tilt angles of the two Gtalons relative to each 
other while the laser output wavelength was moni- 
tored. By this method the Stalon spectral modes were 
aligned to provide for smooth tunability and relatively 
narrow linewidth. The resultant tunable wavelength 
range of the TEMqo, Q-switched Ho laser is shown in 
Fig. 2(b). As can be seen, the tuning was smooth and 
ranged from —2.080 to 2.087 The laser linewidth 
(FWHM) was found to be —0.23 cm - 1 as measured by a 
precision scanning Fabry-Perot interferometer. 

B. DIAL System 

Small portions of the laser beam were sampled by a 
pyroelectric detector and a Ge detector for the moni- 
toring of the laser power and the triggering of the data- 
acquisition system, respectively. The wavelength of 
the Ho laser was monitored with a 0.5-m spectrometer 
used as a spectrograph with the output imaged on a 
2 . 1 -* 4 m sensitive TV video camera (Electrophysics 
Model 7290-06). This arrangement permitted visual 
observation in real time of the laser wavelength over a 
range of —60 cm -1 . 

The output of the laser was collimated to a beam 
diameter of —15 mm by a 5X beam expander. The 
main portion of the laser beam was directed with mir- 
rors and aligned in a coaxial configuration with the 
lidar telescope. The direction of the laser beam was 
steered by a 40-cm X 40-cm mirror. The laser beam 
was transmitted from our lidar laboratory stationed on 
the third floor of the laboratory building and aimed 
horizontally either into the air or toward a building 
—1.1 km away. 

The lidar return signals were collected by a 30-cm- 
diameter Newtonian telescope with a focal length of 
1.75 m and a field of view of —1.2 mrad. The optical 
detector was a 1-mm-diameter, liquid-nitrogen-cooled 


InSb photovoltaic detector (Infrared Associates Model 
IS- 1), which had a cold filter with a 0.35-^m spectral 
bandwidth. The D* of the InSb detector was 3.8 X 
10 11 cm Hz ,/2 /W at 2.1 ftm for a 385-kfl load impedance 
(0.3-MHz bandwidth). For our lidar experiments, the 
detector electrical output was shunted into a 1.3-kfl 
load impedance that raised the detector-amplifier 
bandwidth to 5 MHz but lowered the D* of the InSb 
detector to 2 X 10 10 cm Hz 1/2 /W. The output signal 
from the InSb detector was processed through a 100- 
MHz preamplifier (Analog Module Model 352- 1-S- 
1M) and recorded by a 100-MHz digital oscilloscope 
(LeCroy Model 9400) controlled by a Macintosh II 
computer for signal analysis. 

It may be added that, although we used the liquid- 
nitrogen-cooled InSb detector in all the reported lidar 
experiments in this paper, we also were able to use a 
room-temperature (red-extended) InGaAs photode- 
tector (Epitaxx Model EPX1000GR21) combined with 
a 300-MHz preamplifier (Analog Module Model 313-1- 
S). It was found that the lidar signal-to-noise ratio (S/ 
N) from the room-temperature InGaAs was only ap- 
proximately a factor of 4 less than that from the liquid- 
nitrogen-cooled InSb detector. We believe that ther- 
moelectrically cooling the InGaAs detector and 
reducing the amplifier bandwidth may result in superi- 
or performance by the InGaAs detector over the InSb 
detection; we are currently studying this. 

III. Lidar/DIAL Measurements 

Using the lidar system shown in Fig. 1, we measured 
range-resolved lidar returns from atmospheric aerosol 
backscatter. In addition, we also measured lidar re- 
turns from hard targets to perform path-averaged 
DIAL measurements of atmospheric water vapor. 
The single-pulse S/N of the aerosol backscattered re- 
turn signal was —6-12 at a range of 200 m and varied 
with atmospheric conditions. The single-pulse S/N of 
the return signal from the hard target (brick wall) at a 
range of 1.1 km was —130. 

A. Path- Averaged Water-Vapor DIAL Measurements 

There are several absorption lines of water vapor with- 
in the tuning range of the Ho:YSGG laser that can be 
used for the differential absorption measurement of 
water vapor in the atmosphere. This may be seen in 
Fig. 3, where the computer-generated synthetic trans- 
mission spectrum of the atmosphere within the tuning 
range of the Ho:YSGG laser is shown. 18 The three 
strong absorption lines in Fig. 3 are due to water vapor, 
whereas the other, weaker, lines are due to atmospher- 
ic CO 2 . Because the tuning range of Ho:YSGG is from 
2.080 to 2.089 ^m, the absorption line of water vapor 
near 2.0848 (vacuum wave number 4796.554 cm" 1 ) 
may be used for DIAL measurements. 

The backscattered lidar returns from a brick build- 
ing at a range of 1.1 km were measured as the Ho-laser 
wavelength was tuned across the water-vapor absorp- 
tion line near 4796.5 cm - 1 . Figure 4 shows a plot of the 
normalized lidar returns obtained in 1 June 1990. As 
can be seen, the strong water-vapor absorption line is 
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Fig. 3. Calculated atmospheric transmittance near the 2.084 ^m 
wavelength region over a 10-km horizontal path showing overlap of 
the Ho:YSGG laser-tuning range. The three stronger absorption 
lines are due to water vapor, and the weak background absorption 
lines are due to CO 2 . 



WAVE NUMBER (cm ') 

Fig. 4. Lidar returns backscattered from a building 1.1 km away as 
a function of the Ho: YSGG laser wavelength. The absorption line is 
due to water vapor in the atmosphere. 

easily observed. The water-vapor absorption line 
near 4796.5 cm -1 has an absorption linewidth 
(FWHM)of 0.173 cm -1 and an absorption coefficient 
at line center of 4.70 X 10~Vatm cm at STP. 19,20 The 
measured water-vapor absorption linewidth was esti- 
mated to be 0.25 cm -1 , which was obtained from the 
deconvolution of the Ho-laser linewidth (0.23 cm ') 
and the measured total linewidth (0.34 cm ‘) as shown 
in Fig. 4. It may be noted that the temperature depen- 
dence of the 4795.5 cm -1 water-vapor absorption line 
can be estimated from the published literature, which 
indicates that the coefficient for the temperature de- 
pendence of the air-broadened half-wdith is 0.64. 19,20 
For the experiments reported in this paper, this correc- 
tion term was of the order of 3%, which is small com- 
pared with other uncertainties associated with our ex- 
periments. In addition, the temperature dependence 
of the Boltzmann population factor may be estimated 
from known theory, but previous water-vapor experi- 
ments have shown that direct absorption measure- 
ments are required if one is to quantify this factor 


precisely. 9 Thus the latter factor could not be taken 
into account in our experiments. However, future 
experiments with a multipass absorption cell are being 

planned to measure this value. 

The concentration of water vapor in Fig. 4 was deter- 
mined by using the Beer-Lambert law for attenuation, 

/ = / 0 exp(— 2 aNR), where I is the received power, / 0 is 
the incident power, « is the absorption coefficient, N is 
the concentration, and R is the range. By scaling the 
effective value of a by the ratio of the theoretical to the 
measured linewidths, the data shown in Fig. 4 indicate 
that ~13.2 Torr of water vapor was present in the 
atmosphere. This value may be compared with that 
obtained with a wet-dry bulb thermometer placed out- 
side the window of the laboratory, which indicated 19.0 
Torr of water vapor and an ambient temperature of 
32.8°C. As is evident, the lidar-deduced value is less 
than that actually measured by ~30%. This discrep- 
ancy most probably is due to errors caused by the laser 
linewidth’s being comparable with the water-vapor 
absorption linewidth and, possibly, to variability in the 
spatial distribution of the water vapor. We are cur- 
rently setting up a temperature-controlled, 1000-m- 
long multipass absorption cell to verify the water- 
vapor absorption coefficients and provide laboratory 
calibration of our DIAL measurements. 

B. Range-Resolved Aerosol Backscatter- Water- Vapor 
DIAL Measurements 

The Ho DIAL system was used to measure the range- 
resolved backscatter from naturally occurring aerosols 
in the atmopshere as the laser was tuned on, off-line 
through the water-vapor absorption line near 4796.5 
^m. To demonstrate the utility of the Ho DIAL sys- 
tem, we measured the variation of the atmospheric 
aerosol-water vapor as a function of time. The hori- 
zontal profiles of atmospheric aerosol backscatter and 
water-vapor concentration were obtained by pointing 
the lidar at approximately a +2° elevation angle and 
collecting 100-shot averaged return signals every 10 
min. Each digitized signal was smoothed over a range 
interval of 30 m, and the aerosol backscatter and wa- 
ter-vapor concentration were calculated by using the 
lidar equation 21 and the DIAL equation, respective- 
ly. 22 The measured aerosol backscatter and water- 
vapor concentration horizontal profiles are shown in 
Figs. 5(a) and 5(b), respectively, as a function of time 
during the morning. The aerosol backscatter volume 
and water-vapor concentration are indicated by linear 
gray scales with dynamic ranges of 6 and 4, respective- 
ly. The overlap function and the aerosol attenuation 
were corrected experimentally for the lidar aerosol 
backscatter data by assuming that the aerosol density 
was homogeneous after 9 a.m. 23 It should be noted 
that this correction is a first-order approximation for 
the lidar aerosol backscatter, and a more rigorous tech- 
nique, such as the Klett algorithm, 24 will be assessed 
for our future Ho lidar data. The DIAL measurement 
results, however, are not dependent on this correction 
technique. 

As can be seen from Fig. 5(a), the atmospheric aero- 
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Fig. 5. Ho lidar/DIAL measurements of the temporal variation of (a) the atmospheric aerosol backscattering profile and (b) the atmospheric 
water-vapor concentration profile. The lidar was pointed at approximately a +2“ elevation angle. 



sol backscatter distribution was somewhat variable, 
especially during the morning sunrise (6:40 a.m.) and 
the associated increase in automobile traffic in the 
surrounding area. It should be noted that the range- 
resolved lidar signal shown in Fig. 5(b) is not so dra- 
matic as that of our previous vertical aerosol-profile 
measurements 5 because our present S/N is lower owing 
to the lower power of the tunable Ho laser and a more 
homogeneous distribution of atmospheric aerosol 
backscatter in the horizontal direction. 

The water-vapor data shown in Fig. 5(a) show that 
the water-vapor density was somewhat constant in 
time and homogeneous over the first 200— 700-m lidar 
range. The temperature was measured and varied 
from 24°C at 5 a.m. to 29°C at 10 a.m. The water- 
vapor concentration over the total path length was 
measured to be ~20— 24 Torr. These lidar data essen- 

tially agree with those obtained by a wet-dry bulb 

thermometer measurement. However, we believe 
that in light of our large absolute error, as seen in the 
path-averaged measurement in Fig. 4, such close 


agreement may be coincidental. Further lidar experi- 
ments are being conducted to quantify our results 
better and to reduce these errors. 

IV. Summary and Future Work 

In summary, we have developed a tunable, eye-safe Ho 
lidar-DIAL system, which has been used, for the first 
time to our knowledge, for the preliminary measure- 
ment of the atmospheric aerosol-backscatter and wa- 
ter-vapor concentration profiles. The Ho-lidar aero- 
sol backscatter profile measurements indicate that 
detection of ~1 km can be achieved by using only 5-10- 
mJ/pulse laser pulse energy. The Ho DIAL experi- 
ments indicate that water-vapor concentration profile 
can also be obtained at such a range but with somewhat 
poor (~30%) accuracy. 

Our plans are to evaluate the Ho lidar-DIAL system 
further by increasing the laser power and better con- 
trolling the laser linewidth and spectral purity. Re- 
cently several laser research groups have improved on 
the operation of the Ho laser, and their results have 
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shown (1) fixed single-frequency, narrow-linewidth 
operation, (2) high-power (1-10-J/pulse) operation in 
the non-Q-switched mode, and (3) moderate-power 
(350-mJ/pulse) operation in the Q-switched mode. 
We are currently upgrading our Ho laser to have high- 
er-output power, of the order of 300 mJ, and to be 
injection locked to a single-frequency, tunable Ho os- 
cillator. In this way, we hope to be able to quantify the 
performance of the Ho DIAL system better, especially 
in its spectroscopic utility for water-vapor measure- 
ments. Such quantitative measurements are required 
if one is to be able to compare directly the potential 
usefulness of the Ho DIAL system with current near- 
IR DIAL measurements of atmospheric water va- 
por. 10 ' 25 * 26 

This research was supported in part by the Florida 
High Technology and Industrial Council and by the 
NASA/Langley Research Center. When this research 
was performed, Sungdo Cha was a visiting scientist 
from the Department of Physics, Kangwon National 
University, Korea. The authors are grateful to the 
reviewers for several suggested improvements to this 
paper. 

References 

1. D. K. Killinger and N. Menyuk, “Laser remote sensing of the 
atmosphere,” Science 235, 37-45 (1987). 

2. R. L. Byer, “Diode laser-pumped solid-state lasers,” Science 
239,742-747 (1988). 

3. G. J. Quarles, A. Rosenbaum, C. L. Marquardt, and L. Es- 
terowitz, “High-efficiency 2.09 ^m flashlamp-pumped laser,” 
Appl. Phys. Lett. 55, 1062-1064 (1989). 

4. M. J. Kavaya, S. W. Henderson, E. C. Rusell, R. M. Huffaker, 
and R. G. Frehlich, “Monte Carlo computer simulations of 
ground-based and space- based coherent DIAL water vapor pro- 
filing," Appl. Opt. 28, 840-851 (1989). 

5. N. Sugimoto, N. Sims, K. Chan, and D. K. Killinger, “Eye-safe 
2.1-Mm Ho lidar for measuring atmospheric density profiles,” 
Opt. Lett. 15,302-304 (1990). 

6. P. W. Baker, “Atmospheric water vapor differential absorption 
measurements on vertical paths with a CO- lidar,” Appl. Opt. 22, 
2257-2264 (1983). 

7. R. M. Hardesty, “Coherent DIAL measurement of range-re- 
solved water vapor concentration,” Appl. Opt. 23, 2545-2553 
(1984). 

8. W. B. Grant, J. S. Margolis, A. M. Brothers, and D. M. Tratt, 
“COj DIAL measurements of water vapor,” Appl. Opt. 26, 3033- 
3042 (1987). 

9. E. V. Browell, T. D. Wilkerson, and T. J. Mdlrath, “Water vapor 
differential absorption lidar development and evaluation,” 
Appl. Opt. 18, 3474-3483 (1979). 


10. G. Ehret and W. Renger, “Atmospheric aerosol and humidity 
profiling using an airborne DIAL system in the near IR,” in 
Digest of Topical Meeting on Optical Remote Sensing of the 
Atmosphere (Optical Society of America, Washington, D.C.j 
1990), paper ThA6-l. 

11. R. M. Schotland, in Proceedings of the Fourth Symposium on 
Remote Sensing of the Environment (Environmental Research 
Institute of Michigan, Ann Arbor, Mich., 1966), p. 273-283. 

12. V. V. Zuev, V. E. Zuev, Yu. S. Makushikin, V. N. Marichev, and 
A. A. Mitsel, “Laser sounding of atmospheric humidity: experi- 
ment,” Appl. Opt. 22, 3742-3746 (1983). 

13. N. Menyuk and D. K. Killinger, “Atmospheric remote sensing of 
water vapor, HC1 and CH^ using a continuously tunable 
CO:MgFj laser,” Appl. Opt. 26, 3061-3065 (1987). 

14. L. F. Johnson, J. E. Geusic, and L. G. Van Uitert, “Coherent 

oscillation from Ho :,+ , Yb :{+ and Er 3+ ions in yttrium 

aluminum garnet,” Appl. Phys. Lett. 7, 127-129 (1965). 

15. N. P. Barnes and D. J. Gettemy, “Pulsed Ho:YAG oscillator and 
amplifier,” IEEE J. Quantum Electron. QE-17, 1303-1308 
(1981). 

16. S. W. Henderson and C. P. Hale, “Tunable single-longitudinal- 
mode diode pumped Tm:Ho:YAG laser,” Appl. Opt. 29, 1716- 
1718(1990). 

17. T. Y. Fan, G. Huber, R. L. Byer, and P. Mitzscherlich, “Spec- 
troscopy and diode laser-pumped operation of Tm, Ho:YAG,“ 
IEEE J. Quantum Electron. 24, 924-933 (1988). 

18. R. A. McClatchey and J. E. Selby, “Atmospheric attenuation of 
laser radiation from 0.76 to 31.25 ^m," Environmental Research 
Paper 460 (U.S. Air Force Cambridge Research Laboratories, 
Bedford, Mass., 1974). 

19. L. S. Rothman, R. R. Gamache, A. Goldman, L. R. Brown, R. A. 
Toth, H. M. Pickett, R. L. Poynter, J.-M. Flaud, C. Camy- 
Peyret, A. Barbe, N. Husson, C. P. Rinsland, and M. A. H. 
Smith, “The HITRAN database: 1986 edition,” Appl. Opt. 26, 
4058-4097 (1987). 

20. R. R. Gamache and R. W. Davies, “Theoretical calculations of 
N^-broadened halfwidths of H 2 0 using quantum fourier trans- 
form theory," Appl. Opt. 22, 4013-4019 (1983). 

21. R. T. H. Collis, “Lidar,” Appl. Opt. 9, 1782-1788 (1970). 

22. R. M. Schotland, “Errors in the lidar measurement of atmo- 
spheric gases by differential absorption,” J. Appl. Meteorol. 13, 
71-77 (1974). 

23. Y. Sasano, H. Shimizu, N. Takeuchi, and M. Okuda, “Geometri- 
cal form factor in the laser radar equation: an experimental 
determination,” Appl. Opt. 18, 3908-3910 (1979). 

24. J. D. Klett, “Stable analytical inversion solution for processing 
lidar returns,” Appl. Opt. 20, 211-220 (1981). 

25. E. V, Browell, A. K. Goroch, T. D. Wilderson, and S. Ismail, 
“Airborne DIAL water vapor and aerosol measurements over 
the Gulf stream wall,” presented at the 12th International Laser 
Radar Conference, Aix en Provence, France, 13-17 Aug., 1984. 

26. E. V. Brownell, S. Ismail, and B. E. Grossmann, “Temperature 
sensitivity of differential absorption lidar measurements of wa- 
ter vapor in the 720-nm region,” Appl. Opt. 30, 1517-1524 
(1991). 


20 September 1991 / Vol. 30, No. 27 / APPLIED OPTICS 3943 


Appendix 3 (Reprint of Journal Article) 


a reprint from Applied Optics 

Heterodyne Doppler 1 -Mm lidar measurement of reduced 
effective telescope aperture due to atmospheric turbulence 


Kin Pui Chan, Dennis K. Killinger, and Nobuo Sugimoto 


We performed an experimental study on the effect of atmospheric turbulence on heterodyne and direct 
detection lidar at 1 M m, employing a pulsed Nd:YAG bistatic focused beam l.dar that permitted simultaneous 
heterodyne and direct detection of the same lidar returns. The average carner-to-noise ratio and statistical 
fluctuation level in the lidar return signals were measured in various experimental and atmosphere condi- 
tions. The results showed that atmospheric turbulence could reduce the effective receiverteles^ 
of the 1-pm heterodyne lidar to <5 cm at a relatively short range of -450 m near the ground. The observed ef- 
fective telescope aperture and heterodyne detection efficiency varied during the daj ' " lh * ' f^l 'ous 
turbulence level changed. At this time, we are not able to compare our experimental lidar da a to a rigorous 
turbulence end lid., detection theory which include, independendy ...h. ntd^ 

receiver, and detector geometry. It is interesting to note, however, that the observed '""'Ution of the 
effective receiver aperture was similar In functional form with those predictions based on the heterodyne 
wavefront detection theory by D. L. Fried [Proc. IEEE 55 ,57^67 (1967) and ^ ^‘erodyn, ^e lidar de tectmn 
theory for a fixed monostatic system by S. F. Clifford and S. Wandzura [Appl. Opt. 20, 514-516 ( 1981 )|. W e 
have also applied such an effective receiver aperture limitation to predict the system perform n 
heterodyne Ho lidar operating at 2 irni. 


I. Introduction 

Because of rapid progress in compact and tunable 
solid state lasers in recent years, 1 there has been in- 
creasing interest in the development of all solid-state 
lidar systems. While several demonstration systems 
of direct and heterodyne (coherent) lidar technology 
using solid state lasers at 1-2 urn have been report- 
ed, 26 detailed studies of such lidar systems still lag far 
behind the extensive wealth of information known 
about the more common 10 -mhi C0 2 coherent lidar 
technology. The future development of solid state, 
coherent lidar is particularly significant for remote 
wind velocity and moving target measurements and 
will have a strong impact on satellite based, global 
wind measurements (such as the NASA/Lidar Atmo- 
spheric Wind Sounder) and aircraft windshear sensors 
for the FAA and NASA. 
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The extension of coherent Doppler lidar technology 
from 10- to the new 1-2-mhi wavelength range will not 
only involve the development of new laser, detector, 
and optical systems, but will also involve detailed 
knowledge of the propagation characteristics of the 
lidar beam through the atmosphere. While there have 
been numerous past studies of lidar propagation in the 
turbulent atmosphere, 8-22 a significant number of 
these studies were directed toward 10-Mm C0 2 Doppler 
lidar measurements, with relatively few involved with 
experiments in the 1-2-Mm wavelength range. Al- 
though the theoretical predictions can be scaled to a 1- 
or 2-Mm wavelength range, experimental data are still 
required to validate these theoretical predictions. 

In this paper, we report the first experimental study, 
to our knowledge, of the effect of atmospheric turbu- 
lence on heterodyne and direct detection of lidar signal 
returns at 1 Mm. Measurements were made using a 
short-pulse (8-ns) coherent Nd:YAG bistatic focused 
beam lidar system, which also permitted simultaneous 
heterodyne and direct detection of the same lidar re- 
turn radiation. The average carrier-to-noise ratio 
(CNR) and signal standard deviation were measured 
as a function of the receiver telescope aperture area. 
The effect of atmospheric turbulence was found to 
reduce severely the effective aperture size of the re- 
ceiver telescope for a 1-Mm heterodyne lidar. 
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At this time we are not able to compare our experi- 
mental data to a rigorous atmospheric turbulence and 
lidar detection theory which includes independently 
variable transmitter, receiver, and detector geometry, 
since, to our knowledge, no such general theory exists 
in the published literature. However, it was found 
that the observed limitation of effective receiver aper- 
ture area was similar in functional form to the related 
atmospheric turbulence predictions based on the dis- 
torted wavefront detection theory by Fried 1 - and the 
heterodyne lidar detection theory for a fixed monosta- 
tic system by Clifford and Wandzura. 19 These results 
quantify the effective receiver aperture area for het- 
erodyne detection at 1 and 2 pm and provide a compar- 
ison of the CNRs achieved with heterodyne and direct 
detection technique. It should be emphasized that 
our experimental lidar data were obtained for the lidar 
propagation along a horizontal path, and the laser 
beam was focused on the target. We plan to study the 
more general case of a collimated or infinity-focused 
beam lidar in the immediate future and hope to be able 
to compare these results to future theoretical analysis. 
In addition, the more practical case of a slant-path 
(variable altitude) lidar will also be studied. 

The organization of our paper starts in Sec. II with a 
brief summary of the theories derived by other re- 
searchers for optical heterodyne detection in the pres- 
ence of atmospheric turbulence, to introduce the ap- 
propriate parameters. Section III presents the 
experimental setup and the measured average CNR 
and standard deviation of the fluctuating lidar signal 
returns. In Sec. IV, the experimental results are com- 
pared to the related theories given by Fried and Clif- 
ford/Wandzura and are shown to be similar in func- 
tional form. The parameter values determined from 
the experiments are then used in an analytical treat- 
ment to provide an estimate of the system performance 
for a 1 -Mm Nd and 2-Mm Ho lidar operating in different 
conditions. Finally, Sec. V gives the conclusion. 

II. Related Atmospheric Turbulence and Heterodyne 
Lidar Detection Theory 

In this section we give a brief summary of the perfor- 
mance estimation for heterodyne detection in terms of 
the CNR with an emphasis on the measurement of 
lidar propagation through the turbulent atmosphere. 
This summary is based on previous work by others and 
aims to help explain the terminology used and to out- 
line some of the previous theoretical works which have 
a bearing on the present experimental work. 

In the heterodyne detection system to be discussed, 
it is assumed that the lidar return radiation collected 
by the receiver telescope is collimated, mixed with the 
local oscillator radiation, and focused onto a photodi- 
ode. In the absence of atmospheric turbulence and 
target speckles, it is assumed that the amplitude and 
phase of the signal electrical field, and </>*, are 
uniform on the receiver aperture plane. As such, the 
time-averaged output of the square-law detector at the 
intermediate frequency i-j-, can be given by 


.* = u- -'hv r “ -4-j '• A* 

where Q is the receiving efficiency of a coherent lidar 
system which has been recently shown to be ~0-0.4 for 
a monostatic lidar, 23 rj is the quantum efficiency of the . 
detector which is assumed to be uniform over the de- 
tector area, e is the charge of the electron, hv is the 
photon energy, A r = ird;/4 is the receiver aperture area, 
d r is the receiver telescope diameter, and E\ 0 is the 
amplitude of the local oscillator radiation which, for 
simplicity, is assumed to be uniform on the detector 
plane. In shot-noise dominated heterodyne detection, 
the rms noise current of the square-law detector can be 
approximated by 24 

(i'i) * 2eB(i}eP l /hv) 1 W 

where P, = A r E\J2 is the collected intensity of the local 
oscillator, and B is the electrical bandwidth of the 
detector amplifier. From Eqs. (1) and (2), the CNR 
for heterodyne detection in the absence of atmospheric 
turbulence and target speckles CNRw.o is given by 

CNR, /iU * = An PJ{hvB) t (3) 

where P, « A r E 2 j2 is the collected signal intensity. It 
is seen from Eq. (3) that the CNR for heterodyne 
detection of a nondistorted signal wavefront is propor- 
tional to the receiver aperture area A r . 

In the presence of atmospheric refractive turbu- 
lence, the wavefront of the lidar signal is distorted 
while propagating in the atmosphere. The combined 
effect of atmospheric turbulence on heterodyne lidar 
detection may be described approximately by a pa- 
rameter po called the atmospheric turbulence field co- 
herence length, which can be given by 8,20 

Pt) = ^2.9\k- | W Cl( 2 )( 1 - z//?) 5/:, dzj , (4) 

where k is the wavenumber, C 2 is the refractive index 
structure parameter of the atmosphere, and R is the 
detection range. 

Fried previously developed an analytical model for 
the heterodyne detection of such an atmospherically 
distorted wavefront, 12 in which the extent of the wave- 
front distortion was characterized using only a given 
value of po* His analytical results pointed out that the 
ensemble average of the time-averaged output of the 
square-law detector should be given in the form of 
Eq. (1) but multiplied by a factor F 0 (U) given as 

F 0 (U) = (16/ir ) f u[cos~ 1 u - u( 1 - u 2 ) m ] 
io 

X exp[— (l/2)(l/u) 5/:, ]d“> (5) 

where U = d r /po- Equation (5) is essentially Eq. (5.16) 
in Fried’s paper 12 but has been slightly modified to use 
the definition of po as used in Refs. 8 and 20. The 
average CNR for heterodyne detection in the presence 
of atmospheric turbulence CNRH,av is thus given by 

CNR „,J d nO = (OHO = CNR,, ^00. (6) 

It should be pointed out that Eq. (5) was derived using 
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Fried’s model which includes a first-order estimation 
of optical heterodyne detection and the effect of atmo- 
spheric turbulence in the form of the magnitude of po, 
but with no specific consideration on the propagation 
of the transmitted and backscattered laser light and 
their effects on the determination of the value of po- 
As can be seen in Eq. (6), Fo((/) is essentially a factor 
describing the reduction in the average CNR from the 
effects of atmospheric turbulence. By examining Eqs. 
(5) and (6), we see that for U « 1 (i.e., d r « po), Fo s 1 
and CNR//,av s CNR„,o. In this case, there is little 
effect from atmospheric turbulence. On the other 
hand, for U » 1 (i.e., d r » Po), F„ s (3 poAF) 2 and 
CNRwav a CNR„, 0 (3po/dr) 2 . Noting that CNR„. 0 is 
proportional to A r = ird 2 /4, CNR//, av is readily seen to 
be proportional to (3p 0 ) 2 - This indicates that the CNR 
of heterodyne lidar detection may not depend on the 
aperture size of the receiver telescope but may be 
limited by the atmospheric coherence length p 0 on the 
receiver aperture plane. 

Clifford and Wandzura have calculated CNR reduc- 
tion factor F 0 for a fixed geometry, monostatic hetero- 
dyne lidar but in addition took into account the corre- 
lation effect of the forward and backward lidar 
propagation through the same atmospheric turbu- 
lence. 19 Their results found that for lidar measure- 
ments where the laser beam is focused on the target 

F () 0/> s 11 + (C/3) 2 r, (7) 

where U = d r /po is the same as in Eqs. (5) and (6). It 
may be noted that in the original work of Clifford and 
Wandzura, 19 F 0 was given by F 0 s [1 + (d f /2p n ) 2 J '■ 
We have replaced this with the current form in Eq. (7) 
because the value of po used by Clifford and Wandzura 
is larger by a factor of ~1.5 than that used in this paper. 
It is also seen in Eq. (7) that for U » 1, F 0 ac (3p 0 /d r ) 2 , so 
that CNR h ,»v is limited by the magnitude of p 0 , a 
similar result as that given in Eq. (5). 

To summarize, the above theoretical studies of Fried 
and Clifford/W andzura give physical insight into some 
of the processes occurring in a lidar experiment. 
Fried’s analysis does not aim to predict a value of p 0 (or 
of C 2 ) but states that if a value for p 0 has been mea- 
sured or assumed, its effect on the heterodyne detec- 
tion efficiency is given as in Eq. (5). The analysis of 
Clifford/W andzura calculates the influence of C n on a 
monostatic lidar system and predicts a heterodyne 


detection efficiency reduction factor which is approxi- 
mated by Eq. (7) for a focused transmitter beam. In 
general, however, there are no published theoretical 
analyses on the detection efficiency of a bistatic or 
monostatic heterodyne lidar which allows the practical 
(i.e., experimental) freedom to vary independently the 
size of the transmitter, receiver, and detector. What is 
required for the lidar system design is such a theory 
which would be able to predict the po value on the 
telescope pupil plane of a bistatic or monostatic lidar 
system and its overall effect on the heterodyne detec- 
tion efficiency. Several research groups are now work- 
ing on these theories, and we hope to be able to directly 
compare our experimental work with them when they 
become available. Until this time, however, we believe 
that comparison of our experimental data to the relat- 
ed theories of Fried and Clifford/Wandzura, while not 
exact, will give insight into the physics involved. 

III. Experiment 

Our lidar experiments were carried out employing a 
short pulse, heterodyne Doppler Nd:YAG lidar at 1.06 
pm, which has been described in detail in a separate 
paper. 6 Figure 1 is a schematic of the lidar system. 
The transmitter laser was an injection -seeded, Q- 
switched Nd:YAG laser with a pulse output energy of 
~100 mJ and a pulse length of ~8 ns. The laser beam 
size was measured to be ~6 mm at the exp(— 2) point of 
its maximum intensity. About 98% of the laser pulse 
energy was attenuated before the laser beam was ex- 
panded sixteen times and transmitted to the distant 
target by a 13-cm diam Newtonian telescope. The 
transmitted laser radiation was focused on the distant 
target. Backscattered laser radiation was collected 
using a Cassegrainian telescope with a full aperture 
diameter of 40 cm and an effective focal length of 5 m. 
The two telescopes were put close to each other with a 
separation of ~25 cm between optic axis and aimed at 
the same target. 

In the present experiments, the receiver aperture 
diameter was made variable over a range of 5-20 cm , to 
study the effects of the aperture size on heterodyne 
and direct detection. After collimation with a lens 
whose focal length was 10 cm, the signal beam was 
attenuated with ND filters and split into two beams 
with a 50/50 beam splitter; the diameter of the signal 
beam was ~l-4 mm, corresponding to the receiver 
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Fig. 1. Schematic of short-pulse Nd:YAG lidar 
system at 1 jim providing simultaneous heterodyne 
and direct detection of the same lidar return. 
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telescope diameter of 5-20 cm. One of the signal 
beams was focused and directly detected by an InGaAs 
PIN detector, and the other was mixed with the 
Nd:YAG local oscillator output, focused, and hetero- 
dyne detected by another InGaAs PIN detector. The 
two detectors used for heterodyne and direct detection 
detectors were identical, both with a lOO-^m diameter 
and a rise time shorter than 0.5 ns. The focal length of 
the focusing lens was chosen to be 4 cm for both hetero- 
dyne and direct detection detectors, so that the dif- 
fraction-limited focal spot size (i.e., ~25-100 ^m with 
respect to the 5-20-cm receiver diameter) could be 
equal or less than the detector diameter. The diame- 
ter of the local oscillator output was measured to be ~1 
mm at the exp(-2) point of its maximum intensity. In 
the heterodyne lidar measurements, it was collimated/ 
'expanded by using either a 2X or a 4X beam expander 
to help match the ~l-4-mm diameter of the lidar 
signal beam. 

By frequency tuning the laser diode pumped cw 
Nd: YAG local oscillator, the intermediate frequency of 
the heterodyne signal was variable over a range of 200 
MHz to 1 GHz. The maximum output power of the 
local oscillator was 4 mW. The detector output was 
amplified using a low noise, wide bandwidth preampli- 
fier with a 26-dB gain, passed through a 200-MHz high 
pass filter, full-wave rectified using a microwave detec- 
tor with appropriate capacitor and amplified I00X 
before being digitized with a high-speed, twelve-chan- 
nel, analog-to-digital (A-D) converter in a CAMAC/ 
computer system. Since no square-law detector was 
used, the digitized heterodyne signal was squared prior 
to data processing in the computer to be directly com- 
parable to standard heterodyne detection theory. 24 
The output of the direct detection detector/amplifier 
was digitized with the A-D converter as well. The 
lidar system was operated at a 10-Hz repetition rate, 
and the computerized data acquisition and processing 
were performed in real time. 

The laser beam was transmitted from the lidar sys- 
tem stationed on the third floor of our laboratory 
building and focused onto a brick wall ~450 m away 
and 10 m above the ground. The focused beam size 
was measured to be ~1 cm. Figure 2 shows an example 
of a single pulse, heterodyne detection lidar return as 
displayed/recorded on a fast 1-GHz oscilloscope. As 
can be seen, the lidar signal consists of the 550-MHz 
heterodyne signal superimposed on top of the 8-ns 
long, direct detection envelope signal. What is also 
shown in Fig. 2 is an example of how the receiver 
aperture size affects the heterodyne detection CNR. 
In Fig. 2(a) where receiver aperture diameter d r was 5 
cm, the oscillatory heterodyne signal is much stronger 
than the direct detection (envelope) signal. On the 
other hand, in Fig. 2(b), where d r = 15 cm, almost equal 
heterodyne and direct signal intensity were observed. 

A. Effect of Telescope Aperture Size on CNR 

To measure experimentally the effect of the receiver 
telescope aperture size on the lidar CNR values, the 
average lidar CNR for both heterodyne and direct 
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Fig. 2. Oscilloscope displays showing a single-pulse heterodyne 
signal of lidar return from a hard target at a range of 450 m. The 
receiver telescope aperture size was (a) 5 and (b) 15 cm. Note that 
the 550-MHz heterodyne signal is superimposed on top of the direct 
detection, 8-ns long pulse envelope. 


detection techniques was measured as a function of the 
receiver aperture area. Figure 3(a) shows an example 
of our data obtained when the aperture diameter of the 
aperture stop placed in front of the receiver telescope 
(i.e., A4 in Fig. 1) was changed from 5 to 20 cm. The 
measurement was made at noon on a sunny day in Jan. 
1990 (temperature ~20°C). Each data point shown 
was an average of 600 returned pulses, and the whole 
experiment was carried out in about half an hour. In 
this experiment, the diameter of the collimated local 
oscillator beam was matched to that of the signal by 
varying the diameter of aperture A 2 placed in front of 
the local oscillator. The collection apertures in front 
of the detectors (A1 and A3) were fully open and essen- 
tially not used. 

As seen in Fig. 3(a), the intensity of the direct detec- 
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and direct detection detector. 




tion signal increased in proportion to the telescope 
aperture area due to the effect of aperture speckle 
averaging. 25 However, the heterodyne signal intensity 
is seen in Fig. 3(a) to initially increase in proportion to 
the telescope aperture area but then saturates to a 
maximum value. This maximum value was observed 
to change in time by as much as a factor of -2-3 during 
the day, presumably the result of the changes in the 
level of atmospheric turbulence; this is discussed be- 
low. In the data shown in Fig. 3(a), little improvement 
in CNR was made by increasing the receiver aperture 
size larger than —10 cm. It should be also added that 
despite the complex detection procedure, heterodyne 
detection did not offer superior performance in terms 
of the CNR compared with the straightforward direct 
detection because of the small value of the effective 
receiver aperture size and the wide electrical band- 
width (1 GHz) of our lidar system. 

We also performed a similar experiment but 
changed the diameters of the apertures in front of the 
detectors (A1 and A3), instead of the aperture in front 
of the telescope. This was done to mimic the condi- 
tions assumed for several atmospheric turbulence the- 
ories which take into account a variable aperture stop 
in front of the detector and then backward propagate 
the local oscillator radiation to the target through the 
telescope. Figure 3(b) shows our experimental lidar 
data, which was obtained by varying the diameter of 
the collection apertures placed in front of the detectors 


(A1 and A3) d c over the range of 1-4 mm and leaving 
A2 and A4 fully open. The equivalent receiver tele- 
scope aperture area A r was estimated from A r = n{f r dJ 
2 f { ) 2 where f r and ft are the focal length of the receiver 
telescope and the collimating lens, respectively. As 
can be seen, the results shown are essentially the same 
as those presented in Fig. 3(a). 

During the course of the above experiments, we were 
able to measure an approximate value of the level of 
the atmospheric turbulence using a new instrument 
which acquires the 2-D Fourier transform/power spec- 
trum of the telescope image of a distant light source, 
recorded by using a video camera. 2 ^ Our measure- 
ment indicated that the value of C 2 n varied during the 
course of our lidar experiment from —2 to 7 X 10 
m~ w at 10 m above the ground (horizontal path, tem- 
perature range from 18 to 24°C, sunny day). From the 
approximation of C 2 (h) = C„(l)h 4/ *\ where C;(l) is 
the value at h = 1-m height, 27 we estimate a corre- 
sponding value for C 2 (l) of —4 X 10 -13 to 1.4 X 10 12 
m" 2/3 during the day. We are currently in the process 
of comparing the C 2 n value measured with our video 
camera instrument with that measured with a conven- 
tional scintillometer recently obtained on loan for 
NOAA. In this way, a better quantitative value for C n 
can be obtained. [Note added in Proof: We have 

recently been able to use the NOAA Cn 2 scintillometer 
in our experiments and have found a good agreement 
(within 30%) of the Cn 2 value measured with our video 
camera system.] 

B. Effect of Telescope Aperture Size on Signal 
Fluctuations 

The statistical fluctuation level of the lidar signal 
intensity was measured for both the heterodyne and 
direct detection techniques. Figure 4 shows the nor- 
malized standard deviations of the lidar signal returns 
shown in Fig. 3(a). As seen, the differences between 
the results measured with the two detection tech- 
niques are evident. The normalized signal standard 
deviation (cr) for heterodyne detection did not appear 
to depend on the receiver aperture area in our experi- 
mental conditions and had a value of around 1.1. This 
observed value is slightly larger than unity, which is 
expected for the familiar Rayleigh distribution associ- 
ated with heterodyne detection. 25 On the other hand, 
the value of <j for direct detection measurement is seen 
to decrease as the telescope aperture area is increased. 
This is, as expected, the result of aperture speckle 
averaging. 25 

The value of c in Fig. 4 can be compared to previous 
reports by Holmes et a/. 15 * 18 In their works, the signal 
covariance <x 2 for a nonfocused point receiver (direct 
detection) was found to approach a maximum value of 
about 1.3, when the Rytov logarithm amplitude vari- 
ance a\ = 0A24C 2 n k v *R [{/ * was nearly 0.1. Using typi- 
cal values for our experiment of R = 450 m, A = 1.06 ^m, 
and C 2 n = 7 X 10~ 14 m“ 2/3 , one calculates that a 2 x = 0.05. 
From Fig. 1 of Ref. 18, the predicted value of <x corre- 
sponding to <t 2 = 0.05 is -1.1. This value is in good 
agreement with our heterodyne lidar data shown in 
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Fig. 4. Measured normalized signal standard deviation of hetero- 
dyne (open circles) and direct detection (solid circles) lidar returns 
from a hard target at a range of 450 m as a function of the receiver 
telescope aperture area. 

Fig. 4. In this case, we are comparing our heterodyne 
lidar data to that of a specific direct detection mea- 
surement that employed a nonfocused point receiver, 
because both detection processes are essentially simi- 
lar in detecting an interference (i.e., speckle) pattern. 
On the other hand, our direct detection lidar data 
shown in Fig. 4 also indicate that the value of a may 
approach 1.0 for a point receiver (nearly zero aperture 
area). As can be seen, our values for a are consistent 
with the study of Holmes et al . 

The reduction in signal fluctuation observed for our 
1- M m direct detection lidar in Fig. 4 may be contrasted 
with that normally expected using a 10-pm CO 2 lidar. 
From Eq. (4) it can be calculated that at X = 1.06 pm, po 
is — 1.4 cm, assuming that the atmospheric parameter 
C“ = 7 X 10" 14 m~ 2/:l is uniform along the 450-m long 
lidar propagation length. In the same atmospheric 
conditions, the predicted po at X = 10 pm is —22 cm. 
Because of the smaller atmospheric speckle size at 1 
pm, aperture speckle averaging is more effective for 
direct detection of Nd:YAG lidar returns than at 10 
pm. For example, a normalized signal standard devi- 
ation as low as 0.2 was observed in our 1-pm lidar 
measurement at a detection range of 450 m, while a 
similar value for a was measured with a 10-pm CO *2 
lidar at a much longer (3 km) range. 

C. Temporal Changes of Lidar CNR and Signal 
Fluctuation During the Day 

It is instructive to show the differences in perfor- 
mance between the heterodyne and direct detection 
techniques in different atmospheric conditions. Fig- 
ures 5(a) and (b) depict the average CNR and normal- 
ized signal standard deviation of the hard-target lidar 
returns, recorded between 6:30 a.m. and 12:00 a.m. on 5 
Feb. 1990. The weather condition was sunny, and the 
temperature varied from 18 to 23°C. The detection 
range was —450 m, and the receiver telescope diameter 
was always fixed at 8 cm. Each data point shown is the 
average of 600 lidar returned pulses. In the present 
experimental conditions, the receiver aperture diame- 
ter was several times larger than the value of po, which 
was estimated to vary around 1—3 cm. As seen in Fig. 
5(a), a CNR value of —200 was measured with hetero- 
dyne detection at 6:30 a.m. when the atmosphere was 



TIME 

(b) 


Fig. 5. Measured variation of (a) average CNR and (b) normalized 
signal standard deviation of heterodyne (open circles) and direct 
detection (solid circles) lidar returns from a hard target at a range of 
450 m as a function of time during the day. 

relatively calm. It dropped to —120 at noon when the 
atmosphere became most turbulent. On the other 
hand, the signal intensity measured with direct detec- 
tion was nearly constant during this time period. As 
seen in these experimental results, the heterodyne de- 
tection efficiency is degraded by the presence of atmo- 
spheric turbulence, and this effect changes in time as 
the level of atmospheric turbulence varies during the 
day. 

The result shown in Fig. 5(b) is quite similar to that 
shown in Fig. 4. The normalized signal standard devi- 
ation of the direct detection was measured as —0,8 
when the atmosphere was relatively calm and —0.6 at 
noon when the atmosphere was most turbulent. The 
value for a may be approximated by a = 1/M 172 , where 
M is the ratio of the receiver aperture area to the 
speckle pattern correlation area, 26 which, in the 
present experiment where the laser beam was focused 
on the target, is essentially the number of atmospheri- 
cally induced speckles on the receiver aperture plane. 
As such, the observed differences in <j value are be- 
lieved to be the result of the nearly 30% change in the 
value of po during this time period. 

IV. Analysis and Discussion 

A. Comparison of Experimental Lidar CNR to Related 
Atmospheric Turbulence and Lidar Detection Theory 

At this time, it is not possible to compare our hetero- 
dyne lidar data to a rigorous atmospheric turbulence- 
lidar detection theory which allows the freedom to vary 
independently of the parameters of the transmitter, 
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Fig. 6. Comparison of measured lidar CNR values (open circles) 
and theoretical predictions using the theories presented by Fried 12 
(solid line), and Clifford and Wandzura ,f> (dashed line). 

receiver, and detector diameter of a bistatic lidar sys- 
tem, However, it is shown in this section that our 
experimental data follow closely the functional form of 
the related theories by Fried and Clifford/Wandzura. 

The theoretical average CNR for heterodyne detec- 
tion can be calculated employing Eq. (6), using Eq. (3) 
to obtain CNR^.o* and Eqs. (5) and (7) to obtain F 0 . In 
Eq. (3), the intensity of the lidar signal backscattered 
from the hard target P s may be approximated by 

P, * \PJArK exp(-2 aR)]/R\ (8) 

where P t is the transmitted laser power, f is the target 
reflectivity, A r is the receiver telescope aperture area, 
K is the overall optical efficiency, a is the absorption 
coefficient of the atmosphere, and R is the lidar detec- 
tion range. 

Using the parameter values for our lidar experi- 
ments, the theoretical CNR values were calculated and 
are shown in Fig. 6. In Fig. 6, the solid line was 
calculated using Eq. (5) given by Fried, the dotted line 
used Eq. (7) given by Clifford and Wandzura, and the 
open circles are our experimental data as repeated 
from Fig. 3(a). In these calculations, the estimated 
values for our lidar experimental parameters were: P { 
= 5 mJ/8 ns = 6 X 10 5 W, f = 0.1, K * 0.1 X 0.02 (0.02 
due to 98% signal attenuation with ND filters), 
exp(-2aL) = 0.9, R = 450 m, v = 0.8, hv = 1.86 X 10“ 19 
J,B = 500 MHz, and C 2 = 7 X 10~ 14 m“ 2/3 . For the case 
of heterodyne detection, we estimate that the optical 
detection efficiency is further reduced by a factor of G, 
which is the result of the mismatch in optical align- 
ment (geometrical form, wavefront, polarization, 
etc.) 29 ’ 30 and the heterodyne receiving efficiency fl. 23 
In the experimental data shown in Fig. 6, G was esti- 
mated to be —3.5% by comparing the experimental 
CNR with the theoretical predictions. 

As seen in Fig. 6, the experimental data are similar in 
functional form with the related theoretical predic- 
tions, where po = 1.4 cm was assumed for our experi- 
mental conditions ( C 2 n s7X 10” 14 m” 2/3 and R =s 450 
m). It is also noteworthy that although Eq. (7) was 
derived for a focused beam, monostatic lidar, very 
close agreement was observed with this prediction and 
our experimental data measured with a focused beam, 
bistatic lidar system. This may be partially explained 
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Fig. 7. Calculated results of (a) heterodyne lidar detection CNR in 
the presence of atmospheric turbulence and (b) the normalized value 
to the maximum CNR that a heterodyne lidar achieves with an 
effective receiver telescope aperture area of r = rr(3pn) 74. The 
receiver telescope is varied over a range of 0-700 cm", i.e., telescope 
diameter of 0-30 cm. 

by the fact that Eq. (7) is an approximated form of a 
more rigorous equation and is accurate within 20%. 19 
In addition, in our bistatic lidar system the separation 
of the transmitter and receiver telescope optic axis was 
25 cm, so that there could still be partial correlation of 
the forward and backward lidar propagation. Howev- 
er, a more rigorous analysis as well as further experi- 
mental investigation is required to better quantify 
these results. 

To compare the CNR for direct detection, we as- 
sume that the dominant noise sources are the detector 
dark current and amplifier noises. In this case, the 
voltage CNR for direct detection can be approximated 
by 

CNR/, = P,/(NEP^ + NEP^,) 1 '’, (9) 

where NEPn/„ = (hv/r))(2IpB/e)^ 2 and NEPm = (hv/ 
t})(2 IfrB/e ) l/2 are the noise equivalent powers of the 
detector and amplifier, respectively, Id is the detector 
dark current, and Ie is the amplifier noise current, 24 
Using the estimated values for our detection system of 
NEP DL = 4X 10" 9 W, and NEP^ = 3 X 10“* W, the 
calculated CNRp is found to be —90 for a receiver 
telescope with d r = 8 cm, which is larger than the 
measured value CNR p s 45 by a factor of 2. 

Finally, it should be noted that the fitting of the 
theoretical curves in Fig. 6 to our experimental lidar 
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data was made by adjusting only one parameter, which 
was the value of To illustrate this functional de- 
pendence of CNR on p 0t Fig. 7(a) shows the numerical 
prediction of the heterodyne lidar detection CNR in 
the presence of atmospheric turbulence, using the het- 
erodyne wavefront detection theory by Fried [i.e., Eq. 
(5)]. In Fig. 7, the receiver telescope aperture area was 
varied over a range of 0-700 cm 2 (i.e., d r = 0-30 cm) and 
four values of pu (1, 2, 5, and 10 cm) were chosen to 
represent different levels of atmospheric turbulence. 
As seen in Fig. 7(a), the value of po dramatically 
changes the saturated value of CNR at aperture sizes 
much larger than po. The maximum value for CNR is 
proportional to p 2 , which is explained in the next sec- 
tion. 

Figure 7(b) shows similar data, but the values of 
CNR have been normalized to the maximum value 
CNR max obtained at very large aperture size. As seen 
in these calculated results, the functional curve fitting 
of po to CNR is quite sensitive and appears to yield a 
unique value of po with a small relative error. While 
the calculated maximum CNR can be fitted to the 
experimental result by using the scale factor G, the 
shape of the curve is not changed by the magnitude of 
G but determined solely by po. 

B. Predicted Effective Receiver Telescope Aperture for 
1-pm Coherent Lidar 

Thus far it has been experimentally shown that het- 
erodyne detection efficiency can be severely degraded 
by atmospheric turbulence and is related to the reduc- 
tion in the effective receiver aperture for a coherent 
lidar system. Here effective receiver aperture area 
A e fj is defined to be equivalent to an aperture area on 
which the signal wavefront is uniform, so that Eq. (3) is 
directly applicable for calculating the lidar CNR value. 
It is readily shown from Eqs. (3) and (6) that A t u is 
given by 

A vn = F t} (U)A n (10) 

where A r = 7rd 2 /4 is the receiver telescope aperture 
area. Equation (10) is valid only when the target 
speckle size is much larger than the value of po, so that 
the target speckle effects are negligible. As shown in 
Sec. II, for a receiver telescope whose aperture size is 
far larger than p 0l Fo(U) approaches (3po/d r ) 2 . As a 
result, the maximum value of A e ff will be Aeff.max * 
tt(3po) 2 /4, i.e., a receiver telescope with a maximum 
effective aperture size of d e ff,max = 3po- In our experi- 
ments where po was typically measured to be ~1.4 cm, 
deff t max ^ 4 cm. 

Determining the optimal receiver telescope size for a 
heterodyne Doppler lidar is a very practical question. 
As a design method, the useful diameter of the receiver 
telescope d' r may be determined by the following equa- 
tion: 

CNR #/<UV (*W = (11) 

where 'F is a factor that gives the CNR value of the lidar 
system as a fraction of its maximum CNR value. To 
illustrate the use of Eq. (11), we take three examples of 



Fig. 8. Calculated results of the maximum effective receiver tele- 
scope diameter (solid lines) and the practically useful telescope 
diameter (dashed lines) for a 1.06-/4m Nd:YAG heterodyne lidar as a 
function of detection range. The atmospheric parameters C" - \ X 
10 _,2 ,5 X 10~ M , and 2 X 10~ ,5 m -2/:i are assumed to be uniform along 
the horizontal lidar propagation path at a height of 1, 10, and 100 m, 
respectively. 


>F = 70, 80, and 90%. Substituting Eqs. (3), (5), and (6) 
into Eq. (11) yields the numerical solution of d r s 5po, 
6p 0 , and 9po for 'F = 70, 80, and 90%, respectively. This 
result shows that it may not be practical to use a 
receiver telescope whose diameter is much larger than 
6po, since no significant improvement will be made by 
increasing d r beyond this value. As can be seen, the 
above result is consistent with the experimental data 
and theoretical predictions shown in Figs. 6 and 7. It 
may also be noted that the numerical solution of Eq. 
(11) using Eq. (7) yields nearly the same results. 

The maximum effective receiver telescope diameter 
for a 1.06-pm Nd:YAG heterodyne lidar with d r = 3po 
was calculated using Eq. (4) as a function of the detec- 
tion range, and the results are shown in Fig. 8 (solid 
lines). For references, the calculated results of the 
practically useful telescope diameter d r = 6po, is also 
shown (dashed lines). In Fig. 8, three experimental 
conditions are shown, namely horizontal lidar mea- 
surements at an altitude of 1, 10, and 100 m. In these 
calculations, it is assumed that the atmospheric pa- 
rameter C 2 is uniform along the horizontal lidar propa- 
gation path, with a value of C 2 = 1 X 10 -12 , 5 X 10“ 14 , 
and 2 X 10 -15 m -2/3 at/i = 1, 10, and 100 m, respective- 
ly. As seen in Fig. 8, the maximum effective receiver 
telescope diameter for a ground-based heterodyne 
Nd:YAG lidar at 1.06 pm may be quite small, depend- 
ing on the level of atmospheric turbulence. These 
results also imply that a heterodyne Nd:YAG lidar has 
higher efficiency for airborne measurements. For 
ground-based operation, it will suffer from greater re- 
ceiver coherence loss because of the smaller atmo- 
spheric speckle size. In the latter case, however, one 
may use a 2-D detector array to perform heterodyne 
speckle averaging which will enhance the detection 
CNR and help to overcome this limitation. 31 [Note 
added in Proof: We have also recently used a Monte 
Carlo simulation method 31 to predict the reduced het- 
erodyne detection efficiency due to atmospheric tur- 
bulence, and have found similar results as those shown 
in Fig. 6.] 


2624 


APPLIED OPTICS / Vol. 30. No. 18 / 20 June 1991 




Fig. 9. Calculated results of the maximum effective receiver tele- 
scope diameter (solid lines) and the practically useful telescope 
diameter (dashed lines) for a 2.1-jim Ho heterodyne lidar as a func- 
tion of detection range. The atmospheric parameters are the same 
as those used in Fig. 8. 

C. Predictions for Coherent Doppler 2-pm Ho Lidar 

We have used the above theory (Sec. II) to predict 
the performance of a coherent Doppler 2.1-pm Ho 
lidar. Such a Ho lidar is particularly useful for remote 
sensing because of its eye-safe wavelength and poten- 
tial for diode laser pumping. Figure 9 shows the calcu- 
lated result of d r - 3p 0 (solid lines) and d r = 6p 0 (dashed 
lines) for a heterodyne Ho lidar operated in the same 
experimental and atmospheric conditions as those of 
Fig. 8. By comparing the results in Figs. 8 and 9, it can 
be seen that a larger receiver telescope aperture can be 
effectively employed for heterodyne lidar detection at 
the longer wavelength, since po at X = 2. 1 pm is estimat- 
ed to be ~2.3 times larger than that at X = 1.06 pm. 
For the ground-based heterodyne Ho lidar measure- 
ment at a range of R = 1-2 km, the maximum effective 
telescope diameter d e ff,max = 3po could be larger than 5 
cm. At the higher altitude of h- 100 m, d e rr,ma* * 40 
cm at R = 1 km, and deff.ma* a 10 era at R = 10 km. 
These predictions indicate that the operating perfor- 
mance of a coherent Doppler Ho lidar will only be 
marginally reduced by atmospheric turbulence. Of 
course, further experiments will have to be conducted 
to verify these predictions. 

D. Comparison of Heterodyne and Direct Detection CNR 
at 1 and 1 0 pm 

In our experiments with the short pulse, coherent 
Nd:YAG lidar at 1.06 pm, we found that nearly equal 
CNRs were measured with heterodyne and direct de- 
tection of lidar radiation backscattered from atmo- 
spheric aerosols at a relatively short range of ~500 m. 
To compare this experimental observation to the the- 
ory, we examine the ratio of heterodyne CNR to direct 
detection CNR, T = CNR H , a v/CNR D) which, from Eqs. 
(3), (6), and (9), is given by 

r = (?i/Xi'fl)F 0 (/>i„cf r )G(NEP/ ) /, + NEP 2 ,)" 2 , (12) 

where we have added CNR reduction factor G for 
heterodyne detection, which, as mentioned earlier, is 
to compensate for the signal loss from mismatch in 
optical alignment. Figure 10(a) shows the calculated 
result of T for a ground-based ( h = 10 m, horizontal 
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Fig. 10. Comparison of heterodyne and direct detection ( NR for a 
ground-based (a) 1-^m Nd:YAG lidar and (b) 10-^m C0 2 lidar, both 
using a 20-cm diam receiver telescope. The atmospheric parameter 
C 2 = 5 x 10 -M m" 2 ^ is assumed to be uniform along the lidar 
propagation path and two different electrical bandwidths of 10 and 
500 MHz are shown. 

propagation path) Nd:YAG lidar, where two electrical 
detection bandwidths were taken as parameters, 
namely, fi = 500 MHz for high velocity (1-km/s) Dopp- 
ler measurement, and B = 10 MHz for conventional 
wind sensing. The other typical values were X = 1.06 
pm, C 2 n - 5 X 10 -14 m _2/3 , n = 0.8, Id = 2nA, Ie = 3 pA, 
d r = 20 cm, and G = 0.4. As seen in Fig. 10(a), for B = 
500 MHz, r = 3 at R = 500 m, showing good agreement 
with our experimental result. This also indicates that 
while a wide bandwidth heterodyne lidar at 1 pm is 
feasible for high velocity Doppler measurement, it 
does not provide a much higher detection sensitivity 
compared with that of a direct detection lidar. How- 
ever, as seen in Fig. 10(a), for a conventional Doppler 
wind sensing application where B is ~10 MHz, the 
heterodyne detection technique still offers the advan- 
tage of higher sensitivity at a detection range up to 
several kilometers. 

The result of Fig. 10(a) appears to differ from previ- 
ous experiments of heterodyne detection using 10-pm 
C0 2 laser sources, where the heterodyne detection 
CNR can be several orders of magnitude greater than 
that of direct detection, 24 and range-resolved atmo- 
spheric aerosol lidar measurements using C0 2 hetero- 
dyne lidars have been made at ranges up to 20 km. 32 33 
This is due to the fact that at a 10-pm wavelength, the 
photon energy is ~1/10 smaller and the value of pn is 
approximately sixteen times larger compared with 
those at 1 pm. To illustrate this, Fig. 10(b) shows the 
calculated result of T for a ground-based C0 2 lidar at 
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10 /xm, It is assumed in this example that a liquid 
nitrogen cooled HgCdTe detector is used, with NEP/u 
= 10" 12 W/Hz 1/2 . The other system and experimental 
values are the same as those used for the Nd:YAG 
lidar. It is clearly shown in Fig. 10(b) that the hetero- 
dyne detection technique at the 10-pm wavelength 
offers superior detection sensitivity compared with 
that of the direct detection. The detection sensitivity 
of a C0 2 heterodyne lidar will be only marginally re- 
duced by atmospheric refractive turbulence when the 
detection range is up to 1-2 km. However, for long- 
range remote sensing applications, it is still necessary 
to consider the effects of atmospheric turbulence on 
the effective receiver aperture of a CO 2 heterodyne 
lidar. 

V. Conclusions 

In this paper, we have presented an experimental 
study of heterodyne and direct detection of lidar prop- 
agation through the turbulent atmosphere. It has 
been observed that for heterodyne lidar detection at 1 
jim, serious receiver coherence loss occurred as the 
result of atmosphere turbulence. At a relatively short 
detection range of 450 m, the maximum effective re- 
ceiver telescope diameter was as small as 4 cm. The 
effective receiver telescope diameter was also observed 
to vary during the day as the atmospheric conditions 
and the turbulence level changed. On the other hand, 
as expected, atmospheric turbulence did not affect the 
direct detection CNR. Likewise, because of the small 
atmospheric speckle size at 1 pm, the statistical fluctu- 
ation in the direct detection lidar signal returns was 
significantly reduced by aperture speckle averaging. 
Of course, more extensive experimental studies that 
involve heterodyne lidar measurement at various alti- 
tudes/ranges and direct measurements of C 2 n are re- 
quired to better understand these effects. 

Our observed limitation on the effective receiver size 
for l- pm heterodyne lidar agrees in functional form 
with the related predictions on atmospheric turbu- 
lence and heterodyne detection as given by F ried 12 and 
Clifford and Wanzura. 19 Based on our experimental 
work and these previous theories, we have also made an 
analytical study of the predicted performance of dif- 
ferent lidar systems operated in various experimental 
and atmospheric conditions. It was found that a co- 
herent 1-jxm Nd:YAG lidar may be severely impacted 
by atmospheric turbulence for ground-based measure- 
ments and less so for upper altitude studies. The 
predictions for a coherent 2~pm Ho lidar indicate that 
some reduction is expected in highly turbulent condi- 
tions or when receiver telescopes with very large aper- 
ture size are used. In such cases, some enhancement in 
the heterodyne Doppler measurement performance 
could be gained through the use of a 2-D heterodyne 
detector array to perform speckle averaging at the 
detector plane. 

Finally, it is our hope that our experimental lidar 
data will encourage the development of the appropri- 
ate atmospheric turbulence and heterodyne lidar de- 
tection theory which is directly applicable to a general 


class of lidar systems, so that a more precise compari- 
son of lidar data can be made. 
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Coherent 1-pm lidar measurements of 
atmospheric-turbulence-induced spatial 
decorrelation using a multi-element heterodyne 
detector array 


Kin Pui Chan and Dennis K. Killinger 


Wc have employed a coherent 1-pm Nd:YAG lidar ayalem lo meaaure directly, for the fin.t 
knowledge the reduced spatial coherence length, p„. of the lidar returns caus y 
turbulence. Our experiments were conducted by using a 2 x 2 heterodyne detector array ^jponniUed 
real-time spatial correlation measuremenU of the lidar returns at two d,rrcrent deUc^r spac n^. The 
spatial correlation coefTicienls and spatial «,herence length of the lidar returns from hard Urge were 
measured during a day-to-night time period when the atmospheric turbulence parameter C was 
measured lo vary from 2 x 10 ,s to 2 x 10 " m " 5 . These directly measured values of p, as functio 
C 7 were found to be in good agreement with theoretical predictions. 


I. Introduction 

To date, the majority of coherent Doppler lidars have 
used C0 2 laser sources at 10 pm. 1 In recent years, 
however, there has been an increased interest in the 
development of coherent Doppler lidars that employ 
newly emerging solid-state laser technology, several 
coherent Doppler lidars using 1-2-pm solid-state 
lasers have been developed by different groups. 
Recently we experimentally measured the effect of 
atmospheric turbulence on a coherent 1-p.m lidar 
system and found that the receiver coherence loss 
caused by atmospheric turbulence was much more 
severe at 1-2 p,m than at 10 pm. 6 Such receiver 
coherence loss has been related to the atmospheric- 
turbulence field coherence length, p«, in several opti- 
cal propagation— detection theories. 10 Although the 
characteristics of p„ ore fairly well understood, there 
have been few coherent lidar measurements that 
have measured the value of p 0 directly. Such experi- 
mental information is particularly important for the 
development of these shorter wavelength solid-state 
coherent lidars. 

In this paper, we report, for the first time to our 
knowledge, the direct, real-time measurement of the 
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spatial correlation of the heterodyne (coherent) lidar 
signals under different conditions of moderatc-lo- 
strong atmospheric turbulence. Our lidar measure- 
ments were conducted by using a short-pulse (8-ns) 
coherent Nd:YAG bistatic lidar system and a 2 x 2 
heterodyne detector array, which provided four inde- 
pendent heterodyne lidar signal outputs. The statisti- 
cal cross correlation of the independent signals was 
obtained for the two different detector pair spacings 
and directly yielded the value of p„. This direct 
measurement or p 0 was also compared with the value 
of p„ that was determined from our measurements of 
the atmospheric turbulence parameter C„ . Good 
agreement was observed in these measurements. 

Our spatial correlation measurements are unique 
in that a coherent lidar has been used, but they arc 
similar in functional form to previous lidar or laser 
propagation experiments that used point-source di- 
rect-detection detectors. For example, Pincus et al. 
measured the spatial covariance function or laser 
speckle propagating through the turbulent atmo- 
sphere, using a visible Ar + -ion laser and two spatially 
separated direct-detection detectors. 12 More recently, 
Churnside et al. measured the strong opUcal scintilla- 
tion in the atmosphere, using an Ar* laser and a 
one-dimensional direct-detection detector array. 
What is valuable about the experiments presented in 
this paper is that it is what wc believe is the first l-g,m 
heterodyne lidar experiment that directly measures 
the value of p a under different atmospheric turbu- 
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lencc levels. The values of p 0 were found to be 
consistent with our earlier measurements of the 
effective receiver telescope aperture size for a 1-pm 
coherent Nd:YAG lidar.* 14 

II. Analytical Background 

When lidar radiation is transmitted onto a diffuse 
target, the resultant backscattcrcd speckle lobe will 
be randomly perturbed by atmospheric turbulence ns 
it propagates backward to the transmitter— receiver. 
Although this target speckle-atmospheric turbulence 
interaction is a complex subject, an approximate 
analysis of the relative size of the target speckle and 
the influence of atmospheric turbulence can be made 
as follows. As a first-order approximation, the diame- 
ter of the target speckle, d„ may be calculated by 

d. = 2 \RID, (1) 

where X is the laser wavelength, R is the transmitter- 
receiver to target range, and D is the diameter of the 
transmitted laser beam on the target. As Eq. (1) 
shows, d, is larger at longer distances. In other words, 
the target speckle 6ize becomes larger for longer- 
range coherent lidar measurements in the absence of 
atmospheric turbulence. As a numerical example, for 
a Nd:YAG laser wavelength of X = 1.06 pm and D = 
10 cm, it is calculated that d. = 2 and 20 cm at R = 1 
and 10 km, respectively. 

One of the atmospheric turbulence effects on lidar 
propagation is the reduction of the spatial coherence 
length of the lidar beam. In particular, the atmo- 
spheric turbulence field coherence length p 0 may be 
calculated by 1011 

P# = [ 2.91 k' /; c.'uxi - zinr d7.|r", (2) 


respectively. By comparison of the above-calculated 
results ofr/. and Pn , it is seen that Tor long-range 1-pm 
coherent lidar measurements the strong path-inte- 
grated atmospheric turbulence may be the dominant 
source in the reduction of the spatial coherence 
length of the lidar returns. 

The technique used in this paper i6 to obtain the 
value of p„ directly from the relative spatial cross- 
correlation coefficients of lidar signals that are mea- 
sured with a multielement detector array. Such a 
correlation measurement technique has been used by 
Bincus cl al., 17 who have also theoretically analyzed 
the first- and second-order statistics of laser target 
speckle that is perturbed by atmospheric turbu- 
lence. 1 *' 7 Although their generalized equations cover- 
ing weak-lo-strong path-integrated atmospheric 
turbulence are quite complex, their work also indi- 
cates that for wcak-to-moderate atmospheric turbu- 
lence corresponding to Rytov variance, a x = 
0.124C ~' T/ *R 11 ' 6 . of less than approximately 0.02, 

the cross-correlation coefficient function of the two 
lidar signals, C ; (r), can be approximated for the 
6imple case of a focused transmitted beam as 


C,(r) - (f/,(r,) - (/,)1 f/,(r.) - (/,} l}/"."> 
- cxp|- 2 |pj ~~ 3 ** 

/x 

t Jo 


(4) 


where 7j(r,) and J,(r 2 ) arc the signal intensities that 
are received by the two point detectors that are 
separately located on the receiver aperture plane, r = 
i r - r,| is the detector spacing, and D, is the diameter 
of the laser beam on the transmitter aperture plane. 
From Eqs. (2) and (4) we have 


where k is the wave number and C„ 3 is the refractive- 
index structure parameter of the atmosphere. For the 
case in which C„ 2 is assumed to be uniformly distrib- 
uted along the propagation path, such as in a horizon- 
tal path, Eq. (2) becomes 

p„ = (1.09 A*/? < 3 > 

Equation (3) is often used to deduce the valuc^of p 0 
from the experimentally measured value of C„ ; the 
value of C J is measured most often by using a 
scintillometer. Equations (2) and (3) are valid for 
lidar propagation where the signal transverse coher- 
ence length is larger than the inner scale of the 
turbulence cells, which causes coherence loss. As Eq. 
(3) shows, p 0 is inversely proportional to the detection 
range. In other words, the spatial transverse dimen- 
sions over which the received lidar signal is spatially 
coherent becomes smaller as the detection range is 
increased, so the atmospheric turbulence effect is also 
enhanced. For X = 1.06 pm and a typical value ol 
q J _ 5 x 10 _,< m' ,/:l measured near the ground, p 0 is 
calculated to be ~ 1 and 0.25 cm at R = 1 and 10 km, 


C,(r) = cxp(-2(r/D,)’ - 2(r/p,) ,u |. (5) 

The term or 2 (r/ Po ) s/1 in Eq. (5) differs from that or 
4 ( r / p )*« >n Eq. (47) of Ref. 16 by a factor of 2. This 
difference is, however, merely due to different numer- 
ical definitions of p„ that were derived fromthe 
definitions of the wave structure function U t - 

(r/p„)*'\ used in Rcfs ' 10 and H ’ Qnd D * = 2(r/f> " ) , ’ 
used in Ref. 16, respectively. The value of p„ given by 
Eq (2) is thus larger than that in Ref. 16 by a factor or 
2 VS , i.e., ~ 1.5. Consequently, Eq. (5) agrees with Eq. 
(47) in Rcf. 16. In this study, the value of p„ was 
calculated by using Eq. (2) so that the result would be 
consistent with our previou s r epo rts . 

As can be seen from Eq. (5), if D, is a constant the 
value of p„ can be deduced from the measured values 
of C, at different detector spacings. In Section 111, 
examples of real-time direct measurements of p„ made 
with a 2 x 2 heterodyne detector array are given, 
followed by an analysis of the spatial correlation 
coefficients under different atmospheric turbulence 
levels. 
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Fig. 1. Schematic of l-yim coherent Nd:YAG lidar system employ* 
inga 2 x 2 heterodyne detector array. B.S., beam splitter. 
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III. Experiment 

Experimental Setup 

A schematic of our experimental setup is shown in 
Fig.-l. Our lidar experiments were conducted by 
employing a short-pulse (8-ns) coherent NcfiYAG 
lidar 4 and a 2 x 2 heterodyne detector array. After 
being attenuated to a pulse energy of -5 mJ, the 
output of the injection-seeded, single-frequency, flash- 
lamp-pumped, Q-switched Nd:YAG laser was fully 
expanded and transmitted by a 13-cm-diameter New- 
tonian telescope to a hard target (brick wall) at a 
distance of ~ 450 m. The transmitted laser beam was 
focused onto the target, whose focused beam diame- 
ter was - 1 cm. The lidar return signals were col- 
lected by an 18-cm-diameter ofT-axis parabolic tele- 
scope with a 5-m long effective focal length. The 
transmitter and receiver telescopes were placed in a 
bistatic configuration with a 25-cm separation be- 
tween optic axes. The lidar propagation path was 
nearly horizontal at a height of -6 m above the 
ground. 

Figure 2(a) is a schematic of our heterodyne re- 
ceiver-detector. The lidar signal was collected with 
the receiver telescope and collimated by using a lens 
whose focal length was changed to vary the projected 
detector spacings on the telescope aperture plane. 
The collimated signal beam was then combined with 
the collimated, 1-mm-diameter local-oscillator beam 
by using a 50/50 beam splitter and was not focused 
onto the 2 x 2 detector array. Therefore both the 
lidar signal and the local-oscillator beam were pro- 
jected as a plane-wave (i.e., collimated) beam onto the 
detector surface. Figure 2(b) illustrates the geometri- 
cal dimensions of the detector array, on which four 
identical 1 00- p,m -diameter, 1 -GHz -band-width 
InGaAs photodiode detectors (Epitaxx ETX100) were 
integrated with a 300-fxm cenler-to-center separa- 
tion; the spatial separation of the two diagonal ele- 
ments was 2 x 300 p,m. As Fig. 2 6hows, the projected 
detector spacing on the receiver telescope aperture 
plane, r, may be calculated by using the geometrical 
optics equation r - Tj ; if 0 where r rf is the detector 
spacing on the detector array surface and f r and f c axe 
the focal length of the receiver telescope and the 
collimating lens, respectively. By chosing f c = &-30 
cm in a series of lidar measurements where f r = 5 m 
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Fig. 2. (a) Schematic of the coherent reccivcr-d elector designed 

for the spatial cross-correlation measurement and (b) geometric 
dimensions of the 2x2 detector army. 

and r d = 300 |im, we were able to vary r over the 
range of 0.5-1. 9 cm. 

The heterodyne signal intermediate frequency was 
set at ~ 500 MHz by frequency tuning the laser-diode- 
pumped, cw 4-mW Nd:YAG local oscillator. The het- 
erodyne signal outputs from the four array elements 
were simultaneously acquired in real time by using a 
four-channel signal-acquisition system, in which each 
channel consisted of a 23 dB preamplifier (Analog 
Modules 710-47), a high-pass filter (Mini-Circuit SUP- 
300), a 26-dB amplifier (Mini-Circuit ZFL-500), and a 
microwave detector that was used ns a rectifier 
(Narda 4503). The four-channel signals were trans- 
ferred to a computer-controlled CAMAC system, 
which included a 12-channel variable-gain amplifier 
(Phillips Scientific 7177) and a high-speed 12-channel 
analog-to-digital converter (Lccroy 2249W). Because 
no square-law detector was used, the digitized hetero- 
dyne lidar signal was squared before data processing 
in the computer to be comparable with that of 
standard heterodyne detection theory. 18 The lidar 
system was operated at a 10-IIz repetition rate. 

Results 

Heterodyne lidar measurements were conducted un- 
der a variety of atmospheric conditions. Figure 3 
shows typical results of the heterodyne lidar signal 
intensity that was measured from each of the four 
detector elements of the 2 x 2 array as a function of 
time and that was taken during nighttime. Six hun- 
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Fig. 3. Intensity of heterodyne lidar rcturnB from a hard target as 
a function of time, detected by the individual detector elements of 
the 2 x 2 detector array. 


H 

Z 

UJ 

o 

LL 

LL 

UJ 

O 

O 

z 

o 

h- 

< 

UJ 
0 C 
DC 

o 

0 

1 

CO 

C/) 

o 

DC 

u 


i.o — 

0 . 8 " 

0 . 6 - 

0.4 

0 . 2 ~ 

o.o- 

0 

PROJECTED DETECTOR SPACING ON 
RECEIVER APERTURE PLANE (cm) 


• • 


o 


o 

o 


o 

o 


o 


• C 2 = 2x10 M m 2/1 


O C 2 = 2x10 2/1 


4 


dred lidar return pulses were measured over a time 
period of 60 s. In this measurement, the focal length 
of the collimating lens was chosen as 7.8 cm, so the 
projected detector spacing on the receiver telescope 
aperture plane was estimated to be ~ 1.9 cm and 
-2.7 cm for the 300- and 2 x 300- p.m separated 
detectors, respectively. As Fig. 3 shows, the four 
individual signal outputs from the 2x2 detector 
array, /,-/«, are partially correlated, showing some 
evidence of signal wave-front distortion caused by 
atmospheric turbulence. The average intensities of 
the four channel signals were found to be nearly 
equal, and the statistics of the signals also followed a 
Rayleigh distribution with nearly unity signal stan- 
dard deviation. 15 During this nighttime experiment, 
the atmospheric turbulence parameter C„ was mea- 
sured as C„ 2 = 1.5 x 10- 14 m' 2 ' 3 by a C„ 2 long-path 
scintillometer. 

From the four-channel signals shown in Fig. 3, the 
cross-correlation coefficient C, was measured to be 
0.49 and 0.28 for the projected detector spacing or 1.9 
and 2.7 cm, respectively. These two relative C, values 
can be used in Eq. (5) to obtain a value of p 0 . 
Substituting the experimental data of C, = 0.49 for 
r = 1.9 cm and C,' 3 0.28 for r 3 2.7 cm into Eq. (5) 
yields a value of p 0 s 3.6 cm. 

This experimental value of p 0 can be compared with 
that obtained by using Eq. (3), which relates p„ to C^. 
When the measured values of C„ 2 = 1.5 x 10 m 
and R = 450 m are used, Eq. (3) predicts a value of 
p 0 = 3.6 cm, which is in excellent agreement with that 
measured above with the correlation measurement 
technique that used a 2 x 2 detector array. 

We have also measured the spatial correlation 
values of lidar returns that arc associated with dif- 
ferent atmospheric turbulence levels. Fig. 4 shows 
the measured lidar signal cross-correlation coefficient 
as a function of the projected detector spacing on the 
receiver telescope aperture plane; the projected dctec- 
tor spacing was varied by changing the focal length of 
the collimating lens. Lidar experimental data were 
obtained under two atmospheric-turbulence levels, 
corresponding to C„ 2 = 2 x 1<T 14 m' 2 ' 3 and 2 x 10' 


Fig. 4. Measured croBB-correlation coefficients as functions of the 
projected detector spacing on the receiver telescope aperture plane 
under two different atmospheric-turbulence levels. 

and over a 2-h period when the value of C n 
varied by less than ±10%. As Fig. 4 shows, the 
cross-correlation coefficient C t is, as expected from 
Eqs (4) and (5), much smaller at the higher turbu- 
lence level. For C„ 2 = 2 X KT 14 m’ 2 ' 3 , i.e for a typical 
night-time lidar experiment, the value of C, could be 
reduced by atmospheric turbulence to exp(-l) for a 
relatively small detector separation of -3 cm. For the 
daytime measurement where C„ 2 = 2 x 10 m 
however, a similar value of C, was measured for a 
much smaller detector separation of nearly 1 cm. 

Figure 5 shows a direct comparison of our experi- 
mental result with the analytical predictions for the 
moderate atmospheric turbulence levels correspond- 
ing to C„ 2 s 10’ 14 nT 2 ' 3 . The measured cross- 
correlation coefficients (filled circles) are repeated 



RECEIVER APERTURE PLANE (cm) 

Fig. 5. Comparison of measured cross-correlation coefficients 
(filled circles) and the analytical predictions (solid and dnshed^lincs) 
made by using Eqs. (3) and (4) for d ifTerenl close values of C, . 
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from Fig. 4 /, and the analytical predictions (solid and 
dashed lines) were obtained by using Eqs. (3) and (5) 
to calculate p n and C,(r), respectively; C„ 2 = 1.5 x 
10 -14 , 2.0 x HT 14 , and 2.5 x HT 14 m -2 ' 3 were used. As 
Fig. 5 shows, our experimental data arc in fairly close 
agreement with the analytical predictions. 

It is also found from Fig. 5 that the measured 
results deviate somewhat from the analytical predic- 
tions in such a manner that the measured values of C, 
were slightly smaller than that predicted for the 
detector spacing of less than ~1.5 cm (~p 0 /2) and 
were greater than that predicted for the larger spac- 
ing. This deviation may be due to the finite size 
(nonpoint source) of the 100-p.m-diamcter detector 
elements that were projected onto the receiver plane. 
In our measurements, where f r — 500 cm and f c = 
8-30 cm, the projected finite size of the detector was 
calculated to be in the range of 0.2-0. 8 cm, which is a 
possible contribution to the deviation observed. An- 
other source of measurement error that may contrib- 
ute to this deviation is the variation of C n 2 value 
during the measurement. As can be seen from Fig. 5, 
the functional form of the cross-correlation coeffi- 
cient is determined uniquely by the C n 2 value. Conse- 
quently, a possible ±10 % variation in C„ 2 value 
during our measurements, which typically took ~ 2 h 
to complete, may affect the functional form of the 
measured results to a certain extent. Meanwhile, the 
possible variation in C„ 2 value along the horizontal 
lidar propagation path may also contribute to this 
deviation. 

Finally, the experimental data shown in Fig. 4 
(open circles) for the stronger atmospheric turbu- 
lence level of C„ 2 = 2 x 1(T 13 m" 2/3 were also compared 
with the theoretical prediction of Eq. (4) but showed 
significant deviation. This deviation was not surpris- 
ing, because Eq. (4) is not valid for the strong 
turbulence case where the Rytov variance a x 2 is much 
greater than 0.02. 17 For the experimental data shown 
in Fig. 4 (open circles), <r r 2 was calculated from C„ 
value to be -0.14. Moreover, because the measured 
coherence length of ~ 1 cm was less than the Fresncl- 
zone size (fl\) w2 = 2 cm, the measured result may be 
partially due to intensity fluctuation (scintillation) 
that is not characterized by Eq. (4). Equation (4), 
obtained from Ref. 16, was originally derived to 
include for simplicity only terms involving the phase 
and not the log-amplitude term in the wave-structure 
function. Thus Eq. (4) does not consider the effects 
associated with the Fresnel zone size, (R\) . How- 
ever, these experimental data were found to be simi- 
lar in functional form to the theoretical prediction in 
Fig. 2 of the paper by Holmes el al ., 17 which is valid for 
the strong atmospheric-turbulence levels. 


IV. Summary 

yfc have conducted coherent 1-p.m lidar measure- 
ments of atmospheric-turbulence-induced spatial 
decorrelation under the conditions of moderate-to- 
strong atmospheric turbulence. Our measurement 
method, which employed a 2 x 2 heterodyne detector 


array, was demonstrated to be useful for a simple, 
direct measurement of p n . The spatial cross-correla- 
tion coefficient of the lidar returns from a hard target 
at a range of 450 m was measured to be reduced by 
atmospheric turbulence to a value of exp(— 1) at a 
relatively small projected detector spacing of approxi- 
mately 1 and 3 cm for typical daytime and nighttime 
measurements, respectively. The measured results of 
p 0 were found to be in excellent agreement with the p^ 
value that was deduced from the measurement of C„ 
by a long-path scintillometer . 

Our original intent in using the heterodyne detec- 
tor array was to improve the signal-to-noise ratio of 
the coherent 1-gm lidar system through coherent 
addition of the atmospheric-turbulcnce-distorted li- 
dar returns. 14 As we have 6hown in this paper, an 
added benefit of such a detector array is its ability to 
measure directly the value of p 0 . 
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